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Abstract 

These are the (twice) extended notes of a set of lectures given at the "12th Workshop on Hadronic Inter- 
actions" at the IF/UERJ, Rio de Janeiro (31. 5. - 2. 6. 2000). The lectures aim at introducing essential 
concepts of instanton physics, with emphasis on the role of instantons in generating tunneling amplitudes, 
vacuum structure, and the induced quark interactions associated with the axial anomaly. A few examples 
for the impact of instantons on the physics of hadrons are also mentioned. 



Chapter 1 



Introduction 



Yang-Mills instantons were discovered a quarter of a century ago m . They have furnished the first explicit 
(and still paradigmatic) example of a genuinely nonperturbative gauge field configuration and display a 
wealth of unprecedented geometrical, topological and quantum effects with fundamental impact on both 
ground state and spectrum of nonabelian gauge theories. In fact, one could argue that the discovery of 
instantons marked the beginning of a new era in field theory. In the words of Sidney Coleman : "In the 
last two years there have been astonishing developments in quantum field theory. We have obtained control 
over problems previously believed to be of insuperable difficulty and we have obtained deep and surprising 
(at least to me) insights into the structure of ... quantum chromodynamics." Instantons have remained an 
active, fascinating and highly diverse research area ever since, both in mathematics and physics. 

The following pages contain the notes of a series of elementary lectures on instantons in quantum me- 
chanics and QCD which aim at introducing some of the underlying ideas to a non-expert audience. 



1.1 Motivation and overview 

QCD instantons are mediating intriguing quantum processes which shape the ground state of the strong 
interactions. These localized and inherently nonperturbative processes can be thought of as an ongoing 
"rearrangement" of the vacuum. In addition, there exist many other types of instantons in various areas of 
physics. This became clear shortly after their discovery when it was realized that instantons are associated 
with the semiclassical description of tunneling processes. 

Since tunneling processes abound in quantum theory, so do instantons. We will take advantage of this 
situation by explaining the physical role of instantons in the simplest possible dynamical setting, namely in 
the one-dimensional quantum mechanics of tunneling through a potential barrier. Since the classical aspects 
of the problem are particularly simple and familiar here, we can focus almost exclusively on the crucial 
quantum features and their transparent treatment. Once the underlying physics is understood, this example 
can be more or less directly generalized to the tunneling processes in the QCD vacuum, where we can then 
concentrate on the additional complexities brought in by the field theoretical and topological aspects of 
classical Yang-Mills instantons. We will end our brief tour of instanton physics with a short section on the 
impact of instantons on hadron phenomenology. 

As already mentioned, instantons occur in many guises and we can only give a small glimpse of this vast 
field here. (The SLAC high-energy preprint library lists far beyond 1000 preprints on instantons in the last 
decade.) Moreover, we will spend a large part of the available time on laying a solid conceptual foundation and 
will consequently have less time for applications. To do at least some justice to the diversity of the subject and 
to give an idea of its scope, let us mention a few actively studied topics not touched upon in these lectures: 
instantons have been found in many field theories, ranging from scalar field theory via supersymmetric 
Yang-Mills and gravity to string (or M-) theory. There are interesting relations and interactions between 
instantons and their topological cousins, the non-abelian monopoles and vortices. In several theories, 



probably including QCD, instantons are responsible for spontaneous chiral symmetry breaking. The role 
of instantons in deep inelastic scattering and other hard QCD processes has been examined, and also their 
impact on weak-interaction processes at RHIC, LHC and beyond. In inflationary cosmology and elsewhere 
relatives of instantons (sometimes called "bounces" ) describe the "decay of the false vacuum" . 

Moreover, there are fascinating mathematical developments in which instantons serve as tools. They 
have been instrumental, e.g., in obtaining profound results in geometry and topology, including Donaldson 
theory of four-manifolds (which led to the discovery of infinitely many new differentiable structures on R ). 
Instantons play a particularly important role in supersymmetric gauge, string and brane theories. They 
saturate, for example, the nonperturbative sector of the low-energy effective (Seiberg-Witten) action of 
TV = 2 supersymmetric Yang-Mills theory. 

Theoretical approaches to QCD instantons include well-developed "instanton-liquid" vacuum models p], 
a variety of hadron models at least partially based on instanton- induced interactions (see, e.g. Q), a sum- 
rule approach based on a generalized operator product expansion (IOPE) H, and an increasing amount of 
lattice simulations B. 

There is a vast literature on instantons. In preparing these lectures I have benefitted particularly from 



the articles and books by Coleman H, Kleinert M, Schulman |q], Sakita |J, Polyakov 10 , Vainshtein et 
al [TLlJ, and Zinn- Justin O]. More advanced and technical material on instantons is collected in [113]. A 
pedagogical introduction to instantons in supersymmetric field theories can be found in J14| . 

The program of these lectures is as follows: first, we will discuss the semi-classical approximation (SCA) 
in quantum mechanics in the path integral representation. This will lead us to describe tunneling processes 
in imaginary time where we will encounter the simplest examples of instanton solutions. In the following we 
will turn to instantons in QCD, their topological properties, their role in generating the vacuum structure, 
and their impact on the physics of hadrons. 

1.2 WKB reminder 

As we have already mentioned, instantons mediate quantum-mechanical barrier penetration processes. The 
semiclassical approximation (SCA) is the method of choice for the treatment of such processes. It might 
therefore be useful to recapitulate some basic aspects of this technique, in the form established by Wenzel, 
Kramers and Brioullin (WKB). 

The WKB method generates approximate solutions of the Schrodinger equation for wavefunctions with 

typical wavelength A small in comparison to the spacial variations of the potential. This situation corresponds 

to the semiclassical limit h — ► where 

2irh , , 

A = >0. (1.1) 

P 
Macroscopic structures, for example, behave normally according to classical mechanics because their wave- 
functions vary extremely rapidly compared to the variations of the underlying potentials or any other length 
scale in the problem, even if of microscopic (e.g. atomic) origin. The semiclassical limit is therefore the 
analog of the geometrical (ray) limit of wave optics, where the scatterers are structureless at the scale of the 
typical wavelengths of light and where wave phenomena like diffraction disappear. 

In order to set up the semiclassical approximation, the WKB method writes the wave function as 

$ (a) = e <*(»)/\ (1.2) 

Since $eC this ansatz is fully general. In scattering applications <& is often called the "Eikonal" . Inserting 

+ V (x) - E 



(1.2) into the Schrodinger equation 

2m dx 2 
one obtains the WKB "master equation" 



ip(x) = (1.3) 



*-#*(*)-*.#> M = E-V&. (1.4) 



This nonlinear differential equation can be solved iteratively. One expands 

$ (x) = $o 0) + H®i (x) + h 2 <S> 2 (x) + ..., 



(1.5) 



(the series is generally asymptotic) and equates terms of the same order of H on both sides. For each power 
of H one thereby obtains a corresponding WKB equation. The zeroth-order equation has the solution 



$o (x) = =b / dx'p {x ) with 



p (x) = y/2m [E~V (x)] 



(1.6) 



where p (x) is the momentum of the particle in a constant potential U x — const whose value equals that of 
V at x. If V (x) varies slowly compared to i/> (x), ip indeed experiences a locally constant U x and the O (?i ) 
solution 

V>0 (x) = e */r>Fdx' v /2m[E-V(x>)] ^^ 

becomes a useful approximation. By inserting it back into the Schrodinger equation ( |l.3| ), 

itl V (x) ip (x) 



w +v w- JS ' 



^o (x) = - 



2 ^2m [£ - V (x)] : 



(1.8) 



we confirm that it is a solution up to the correction term on the right-hand side, which is indeed of O (H) 
and proportional to the variation V of the potential. 

The quality of the zero th-order approximation can be gauged more systematically by checking whether 
the neglected term in (|l.4| ) is small, i.e whether 



ih fall 
2m 



J_$'2 

2m 



l cjy2 



d h 
dx $' 



<1. 



(1.9) 



With A (x) = 2-kH/p (x) = 27tS/$q and with <!>(, ~ $' (to leading order) this turns into 



s- <1 



1 


dA 


1 




1 


5A 


2tt 


dx 


2vr 


A 


" 2tt 


A 



(1.10) 



where 5X is the change of A over the distance of a wavelength. From the above inequality we learn that the 
semiclassical expansion is applicable in spacial regions where the de-Broglie wavelength 1) is small compared 
to the typical variations of the potential (therefore highly excited states behave increasingly classical) and 
where it 2) changes little over distances of the order of the wavelength. (The latter condition alone is not 
sufficient since we have compared two terms in a differential equation, not in the solution itself.) 

The semiclassical treatment of tunneling processes, which will be a recurrent theme throughout these 
lectures, involves a characteristic additional step since tunneling occurs in potentials with classically forbidden 
regions, i.e. regions x € [x/,x u ] where E < V (x) . The boundaries x;,x u are th e cla ssical turning points. 
Inside the classically forbidden regions, p (x) becomes imaginary and the solution (1.7) decays exponentially 

H3 ' , / 

funnel (x) = e^™? « >&&&=* . (1 . n) 



A comparison of the two solutions (1.7) and (1.11) shows that, formally, the meaning of allowed and forbidden 
regions can be interchanged by the replacement 



t 



E^iE. 



V^iV 



(1.12) 



which indeed transforms (1.7) into the tunneling amplitude (1.11). In other words, we can calculate tunneling 
amplitudes in SCA by the standard WKB methods, analytically continued to imaginary time. Below we will 
see that this works the same way in the path-integral approach. 



1 Since Gamov's classic treatment of nuclear a-decay as a tunneling process through the Coulomb barrier of the nucleus this 
amplitude is also known as the "Gamov factor" . 



Although the WKB procedure is simple and intuitive in principle, it becomes increasingly involved in more 
than one dimension (except for special cases |l5j) and at higher orders. The O (h n ), n > corrections are 
calculated by matching the solutions in the classically allowed and forbidden regions at the classical turning 
points. This leads to Bohr-Sommerfeld quantization conditions for p (x) when integrated over a period of 
oscillation in the allowed regions. However, the generalization of this procedure to many-body problems and 
especially to field theory is often impractical |16[. Fortunately, there is an alternative approach to the SCA, 
based on the path integral, which can be directly applied to field theories. With our later application to 
QCD in mind, we will now consider this approach in more detail. 



Chapter 2 



Instantons in quantum mechanics 



2.1 SCA via path integrals 

One of the major advantages of the path integral representation of quantum mechanics is that it most 
transparently embodies the transition to classical mechanics, i.e. the semiclassical limit. In the following 
chapter we will discuss this limit and the associated SCA in the simplest dynamical setting in which instantons 
play a role, namely in one-dimensional potential problems with one degree of freedom (a spinlcss point 
particle, say). 

2.1.1 The propagator at real times: path integral and SCA 

The key object in quantum mechanics, which contains all the physical information about the system under 
consideration (spectrum, wavefunctions etc.), is the matrix element of the time evolution operator 



x ! 



e -iHT/h 



(2.1) 



i.e. the probability amplitude for the particle to propagate from 



T T 

Xi at t = — — to xf at t = — . (2.2) 

(In quantum field theory this matrix element, which is sometimes called the quantum- mechanical prop- 
agator, generalizes to the generating functional.) Its path integral representation is 

(x f \e- iHT ' h \ Xl )=N JD^e^ =M [ D[x] {x{ _ T/2)=Xtlx{T/2)=Xf} ei{-T% M ^ (2.3) 



where 

r T/2 r T/2 cm -i 

S[x]= dtC(x,x)= dt \ — x 2 (t) - V [x (t)] \ (2.4) 

J—T/2 J—T/2 ^2 > 

is the classical action along a given path. The integration over the paths can be defined, for example, by 
discrctizing the time coordinate (Trotter formula) into intervals At so that t„ = nAt and 

D[x] := lim — — TT — — dx n 2.5 

L J N~*oo\2TTiht) J- 1 \2nitUj v ' 

where x n = x(t n ) . An alternative representation of D[x] in terms of a complete set of functions will be 
introduced below. 




Figure 2.1: The two paths from {x, t) — (0, 0) to (1, 1) whose action is compared. 



In the semiclassical limit, i.e. for h — > 0, the classical action can become much larger than h. As a 
consequence, the path integral is dominated by the paths in the vicinity of the stationary point (s) x c i (t) of 
the action (if they exist), which satisfy 



_S_ 
5x{t) 



S \x\ = 



dV 

dx 



= 



with the boundary conditions 



T 

■'■/I "2 



T 

■'-' ' 2 



Xf. 



(2.6) 



(2.7) 



These classical paths are important not because they themselves give dominant contributions to the 
integral (in fact, their contribution vanishes since the set of classical paths is of measure zero) but rather 
because the action of the neighboring paths varies the least around them. Thus, infinitely many neighboring 
paths, which lie in the so-called "coherence region" with similar phase factors exp(iS/H), add coherently 
For the paths outside of the coherence region the phases vary so rapidly that contributions from neighboring 
paths interfere destructively and become irrelevant to the path integral. This exhibits in a transparent and 
intuitive way the relevance of the classical paths in quantum mechanics. 

To get a more quantitative idea of the coherence region, let us approximately define it as the set of 
paths whose phases differ by less than 7r from the phase of the classical path (the "stationary phase"), which 
implies 

SS [x] = S [x] - S [xci] < irh. (2.8) 

Thus, for a macroscopic particle with 



S ~ 1 erg sec 



gem" 

I ^-sec 

sec^ 



io 27 h 



(2.9) 



only an extremely small neighborhood of the classical path contributes, since 5<f> = SS/h gets very sensitive to 
variations of the path. A numerical example (borrowed from Shankar J17J) makes this more explicit: compare 
two alternative paths from [x, t) = (0,0) to (x,t) = (1cm, Is) (see Fig. |2.l|) for a noninteracting particle: 
the classical path x c i (t) — t with S [x c i] = J dt(m/2)v [x c i] — mem 2 / (2 sec 2 ) and the alternative path 
x a it (t) — i 2 with S [x a it] = 2m cm 2 / (3 sec 2 ). For a classical particle with 



m = Ig 



5S = S [xait] - S [x c i] 



6 sec 2 



1.6 x 10 26 h 



1.6 x 10 2t> rad ;» n 



the alternative path is extremely incoherent and therefore irrelevant, while for an electron with 



m = 10~ 27 .g 



5S~U 
6 



- rad <C 7r 
6 



(2.10) 



(2.11) 




Figure 2.2: A typical tunneling potential with nondegenerate minima. The total energy (horizontal line) is 
smaller than the hump so that a classically forbidden region exists. 



it is well within the coherence region and makes an important contribution to the path integral! 

Obviously, then, the movement of a free electron (even a very fast one) cannot be described classically. 
One has to resort to quantum mechanics where the electron's path is much more uncertain and, in fact, not 
an observable. There also exist microscopic situations, on the other hand, where the quantum fluctuations 
do not totally wash out the classical results and which can therefore be treated semiclassically. Typical 
examples are the scattering off a slowly varying potential (Eikonal approximation) or the highly excited 
electronic orbits in atoms (Rydberg atoms). 

However - and this is a crucial point for our subsequent discussion - the stationary-phase approximation 
fails to describe tunneling processes! The reason is that such processes are characterized by potentials with 
a classically forbidden region (or barrier) which cannot be transgressed by classical particles. In other words, 
there exist no classical solutions of (2.6) with boundary conditions corresponding to barrier penetration in 
a potential of the type shown in Fig. 2.2, As a consequence, the action S [x] has no extrema with tunneling 
boundary conditions, and therefore the path integral (2.3) has no stationary points. 

Does this mean that the SCA becomes unfeasible for tunneling processes? Fortunately not, of course, as 
the standard WKB approximation shows. The problem is particular to the stationary-phase approximation 
in the path integral framework, which indeed ceases to exist. The most direct way to overcome this problem 
would be to generalize the familiar saddle-point approximation for ordinary integrals. This method works by 
deforming the integration path into the complex plane in such a way that it passes through the saddle points 
of the exponent in the integrand. And since there indeed exist complex solutions of the combined equations 
(2.6), (2.7) in tunneling potentials (in other words, the complex action still obeys the Hamilton- Jacobi 
equation) , a generalization of the complex saddle-point method to integrals over complex paths would seem 
natural. But unfortunately integration over complex paths has not been sufficiently developed (although 
a combined approach to both diffraction and ray optics in this spirit was initiated by Keller |lq| , and 
McLaughlin |H| has performed a complex saddle point evaluation of the Fourier-transformed path integral. 
See also the recent work of Weingarten |2(J on complex path integrals.). 

However, in many instances - including our tunneling problem - an analytical continuation of the path 
integral to imaginary time, as already encountered in the WKB framework in ( 1.12J ) , serves the same purposes. 
This is the way in which we will perform the SCA to tunneling problems in the following sections. 



2.1.2 Propagator in imaginary time 



As just mentioned, the paths which form the basis of the SCA to tunneling processes in the path-integral 
framework can be identified by analytical continuation to imaginary time. This approach was apparently 



first used by Langer (22| in his study of bubble formation processes at first-order phase transitions^. In the 
context of the SCA, tunneling problems in field theory, and instantons it was introduced by Polyakov et al. 
[g, g4[, building on a suggestions by Gribov. 

In order to prepare for the use of this approach, let us first see how the ground-state energy and wave- 
function in tunneling situations can be obtained from the matrix element 

Z{x f ,Xi) = (x f e~ HTE ' h Xi ) (2.12) 



which is obtained from (2.1) by analytically continuing 

t->-ir (=*• T -» -iT E ) . (2.13) 

This procedure (together with its counterpart p — » ip^ in momentum space) is often called a "Wick rotation" . 
(Note, incidentally, that imaginary- "time" evolution as generated by e~ HTB ^ h is not unitary and therefore 



does not conserve probability. It might also be worthwhile to recall that the matrix element ( 2.12 ) plays 
a fundamental role in statistical mechanics (its trace over state space, corresponding to the sum over all 
periodic paths, is the partition function) although this is not the angle from which we will look at it in the 
following.) 

The matrix element Z (xf, Xi) has the energy representation (in the following we drop the subscript E of 
T E ) 

z(x f ,xi)=j2 e ~ EnT/h ^f^ W 3 *) ( 2 - 14 ) 

a 

in terms of the spectrum 

H\n)=E n \n), 1 = ^ |n) (n| (2.15) 

n 

of the (static) Hamiltonian H. The energies E n are the usual real eigenvalues, which remain unaffected by 
the analytical continuation. For large T the ground state dominates, 

Z (x f , x t ) -> e- E ° T ' h (s/10) (Ola:,) , (2.16) 

and the ground state energy becomes 

E = -h lim -lnZ{x f ,Xi). (2.17) 

In order to calculate the ground state energy (and wave function) of a quant u m- mcc hnical system we therefore 
just need to take the T — ► oo limit of the imaginary-time matrix element ( 2.12J ). This convenient way to 



obtain ground-state properties is used, e.g., in Euclidean lattice QCD to calculate hadron masses. In the 
following section we will show how to calculate the matrix element Z(xf,xi) semiclassically in the path 
integral framework. 

2.1.3 Path integral in imaginary time: tunneling in SCA 



The quantum mec han ical propagator (2.12) in imaginary time has a path integral representation which can 



be obtained from (2.3) by analytical continuation: 

Z{x f ,Xi)=N J D[x}e-^ =N [D[x] M _ T/2) = Xi]x{T/2) = Xf} e-*!-T% dTCE ( x > A K (2.18) 



(Note, by the way, that this path integral is much better behaved than its real-time counterpart (and thus 
can be formalized as an integral over generalized stochastic Wiener processes) since the oscillating integrand 
is replaced by an exponentially decaying one.) 

1 Another early application of similar ideas is due to W.A. Miller ]23| . (I thank Takeshi Kodama for drawing my attention 
to this reference.) 



For future applications it is useful to define the "measure" of this integral more explicitly. To this end, 
we expand x (r) into a complete, orthonormal set of real functions x n (r) around a fixed path x (r) as 



x (t) = x (t) + 77 (t) 



and 



77 2 



-T/2 



drx n (r) x m (r) = <S„ 



where ?/ (r) = ^J c c i„ (r) 

^ in (r) x„ (r') = <5 (r - t') . 



Furthermore, a; is assumed to satisfy the boundary conditions implicit in the path integral, i.e. 

x (±T/2) = Xf/Xi, X„ (±T/2) = 0. 

(Note that (for now) there are no further requirements oni.) As a consequence we have 

dc„ 



D[x]=D [ V ] = [] ■ 



'2irh 



(2.19) 
(2.20) 
(2.21) 
(2.22) 



where the normalization factor is chosen for later convenience. 



Now let us get the explicit form of Ce (x, x) (the subscript E stands for "Euclidean time" , which is used 
synonymously with "imaginary time" in field theory) in (2.18) by analytical continuation of the integration 
path as 



iS = i 
After substituting t — 



T/2 



dt 



V\x 



T/2 



/m 



according to (2.13|) we obtain 



dt 



,772 
i (— i) / dr 

J -T/2 



(■ 



?(f'-'"' 



T/2 



= —Sf 



V[x\) =- J dTC E [x] 

' J-TI2 



(2.23) 



(2.24) 



(note that in Eq. ( 2.24 ) and in the following the dot indicates differentiation with respect to r) from which 
we read off 



C E [x] 



-x 2 + V (x) 



(2.25) 



Thus, besides turning the overall factor of i in the exponent into a minus sign (as anticipated in ( |2.18 )), the 
analytical continuation has the crucial implication that the potential changes its sign, i.e. that it is turned 
upside down! 

The consequence of this for the SCA is immediate: now there exist (one or more) solutions to the 
imaginary-time equation of motion 

j- n rS E [x]=mx c t-V'(xc l )=0 (2.26) 

ox (r) 

with tunneling boundary conditions. These solutions correspond to a particle which starts at Xi, rolls down 
the hill and through the local minimum of — V (x) (which corresponds to the peak of the classically forbidden 
barrier of +V (x)) and climbs up to cc/, according to the boundary conditions 



T 

x I ~ 77 I — x ii 



T 
' ' 2 



x f . 



(2.27) 



For later reference we note that the solutions of (2.26) carry a conserved quantum number, the "Euclidean 
energy" 

E E = ^x 2 -V{x) (2.28) 



which follows immediately from Ee = x c i [mx c i — V' (x c i)] and the equation of motion (2.26). 

The essential property of these solutions in our context is that they are saddle points of the imaginary- 
time path integral ( |2.18 ). For h — > the only nonvanishing contributions to the path integral therefore come 
from a neighborhood of x c i (since those are least suppressed by the Boltzmann weight exp (— Se/K)) and 
can be calculated in saddle-point approximation. 
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2.1.4 Saddle point approximation 

Now let us perform the saddle-point approximation explicitly. To this end we write 

x (r) = x c i (t) + 77 (r) 



(2.29) 



and expand the action around the stationary (resp. saddle) point x c i (later we will sum over the contributions 
from all saddle points) to second order in 77, 



S E [x] = S E [xd] + x dl ~ rfr ' 7 ? ( T ) 



S 2 S E [x c 
Sx (t) Sx (r') 



-r,{ T ') + 0{rf) 



= S E M + \ j drr, (r) F (x cl ) 17 (r) + O {r, 3 ) , 

1 J-T/2 

(the first derivative vanishes since Se is minimal at x c i) where we have abbreviated the operator 



S 2 S E [x] 
Sx (r) Sx (r') 



dr 2 



d 2 V(x cl ) 
dx 2 



S{t-t') = F(x c1 )S{t-t') 



(2.30) 
(2-31) 

(2.32) 



which governs the dynamics of the fluctuations around x c i . 

We now expand the fluctuations 77 (r) into the (real) eigenfunctions of F, 



with 



■q (t) = ^2 c nX n (r) , (2.33) 

n 

F (x c i) x n (r) = X n x n (t) . (2.34) 

The x n satisfy the boundary conditions (2.21) and are normalized according to ( 2.20J ) . (Since F is real- 
hermitean and bounded, it has a complete spectrum.) For the moment we will also assume that all eigenvalues 
are positive, A„ > 0. (This is not true in general. Below we will encounter examples of vanishing eigenvalues. 
In the case of so-called "bounce" solutions even negative eigenvalues do occur J2|.) The action can then be 
written as 

S E [x] = S E M + ^ A « c » + ° fa 3 ) ■ ( 2 - 35 ) 



Now we use the definition ( [2.22 ) of the measure to rewrite 



Z(x f ,Xi)=M I D[x]e-^ ~ Ne-^^ / D[ V ]e-^ J-%* d ^(r)H-cMr) 



as 



s fi 1 _— _ r°° 
Zixj^^Ne-^^W 

„ J — OO 



dc 



n 



-AE„a„ 



= Me 



n 



x, v 2nh 

e -2fe A " c -, 



OO 



dc 



II 



'2-Kh 



(2.36) 



(2.37) 



(2.38) 



and obtain, after performing the Gaussian integrations (which decouple and can thus be done independently) , 

(2.39) 



Z(x f ,x l )=JVe-^Y[\- 1 / 2 



AAe-^fr^det^M) 



-1/2 



(2.40) 
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where a sum J^ is implied if there exists more than one saddle point. The formula ( 2.4C ) encapsulates the 
SCA to Z (xf,Xi) up to O (h). 

Above, we have introduced the determinant of a differential operator as the product of its eigenvalues 

detO-J|A„, for 6^ n (x) = \ n i/) n (x) , (2.41) 



which generalizes the standard definition for quadratic matrices. In Subsection 2.3 we will give a more 
detailed definition of such determinants (which renders the above (infinite) product finite by adopting a 
specific choice for the normalization factor Af ) and show how they can be calculated explicitly. 

Let us summarize what we have accomplished so far. The seemingly artificial analytical continuation to 
imaginary times has allowed us to identify those paths whose neighborhoods give the dominant contributions 
to the path integral for a tunneling process in the scmiclassical limit, and to evaluate this path integral 
to O (h) in the saddle-point approximation. In more physical terms the situation can be described as 
follows: for tunneling problems there exist no minimal-action trajectories (i.e. classical solutions) with the 
appropriate boundary conditions in real time. Therefore all trajectories between those boundary conditions 
(over which we sum in the real-time path integral) interfere highly destructively. Still, their net effect can be 
approximately gathered in a finite number of regions in function space, namely those in the neighborhood 
of the saddle points in imaginary time. In other words, while the tunneling amplitudes would have to be 
recovered at real times from a complex mixture of non-stationary paths (a forbidding task in practice) , they 
are concentrated around the classical paths in imaginary time, and are therefore accessible to the saddle- 
point approximation. The destructive interference at real times leaves a conspicuous trace, however, namely 



the exponential suppression of (2.40) due to the Gamov factor exp(— Se/K), which is typical for tunneling 
amplitudes. 



2.2 Double well (hill) potential and instantons 

In the following sections we will apply the machinery developed above to simple tunneling problems with 
potentials which resemble as closely as possible the situation to be encountered later in QCD. 



2.2.1 Instanton solution 

Let us therefore specialize to tunneling processes between degenerate potential minima. A simple potential 
of the appropriate form is 

V(x) = ^(x 2 -x 2 ) 2 (2.42) 



(see Fig. 2.3, the coefficient is chosen for later convenience) which has three saddle points (with finite action, 



see the comment at the end of this section). Two of them are trivial 

1) x c i (r) = xo (2.43) 

2) x c i (r) = -x (2.44) 

and do not contribute to tunneling (since they cannot satisfy the corresponding boundary conditions; they do, 
however, contribute to Z (xq, xq) or Z (— Xo, — Xo))- The third saddle point is time-dependent and corresponds 

12 




Figure 2.3: The double-well potential with %-Q = 10 and xq = 1. 



to the tunneling solution which interpolates between both maximag of — V, 

x cl (-T/2) = +x 0: (2.46) 

x cl (T/2) = -so, (2.47) 

or in the opposite direction, i.e. starting at — xq and ending at xq. 

The above classification of solutions is common to all potentials which look qualitatively like (2.42), and 
all the qualitative results obtained below will apply to such potentials, too. We have specialized to (2.42) 
only because for this choice the tunneling solution can be obtained analytically. The easiest way to get this 
solution starts from the conserved Euclidean energy (2.28) with Ee = 0, which corresponds to the limit 
T — > oo we are interested in (since then the particle in the mechanical analog system starts with (almost) 
zero velocity at xq and thus will need (almost) infinite time to reach — xq): 



E k = —x 2 



V = 



T 



(as usual, the positive square root is implied) or, after separation of variables, 

dx 



T 



2 vn^) 



dr, 



which can be integrated to become 

r x " l ^ T ) dx -*•„ r x ci( T ) 



=F 



_x 



xci(ro) y/V{x) 



dx 



Xct(ro) ^0 



=p arctan h 

a xq 



Xcl(r) 

Xcl(To) 



dr 



(2.48) 



(2.49) 



r . (2.50) 



We now choose to, the integration constant, to be the "center" of the solution in imaginary time by requiring 
x c i (to) = 0, and finally obtain 



3) x c i (t) = xjfj (t) = =fx tanh a(r - r ) . 



(2.51) 



The solution with the minus sign is called the instanton (see Fig. 2.4) of the potential ( |2.42 ). Tunneling 
from — xq to lo corresponds to the anti- instanton solution Xj{t) — —xi (r). Obviously, this nomenclature 

2 Note that, in the exact sense, this solution exists only for T — ► co (which is all we need to extract the ground-state 
properties). At large but finite T the corresponding solution requires a slight change in the boundary conditions to 

x cl (±T/2) = Txo±e (T) , (2.45) 

where e is small. 
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Figure 2.4: The instanton solution in the double well potential of Fig. 2.3 



is mere co nven tion. (We repeat that, from the conceptual point of view, there is nothing special about the 
potential ( 2.42J ). Similar potentials with degenerate minima will have similar instanton solutions although 
those generally cannot be found analytically.) 

Note that instantons are by necessity time-dependent since they have to interpolate between different 
minima of the potential. The name "instanton" (which goes back to 't Hooft) indicates furthermore that 
the tunneling transition happens fast, i.e. almost instantaneousn. Indeed, from the equation of motion (for 
Ee = 0) for the instanton, 

dxi 



-V(xi), 



(2.52) 



we obtain at large r (where V (xi) can be expanded around — xq with V (— Xo) — V (—Xo) = and 
V" (-x ) = 4a 2 m) 



— ~ -2a[x! - (-x )}- 



Separating variables, 



dxi 



d In (xi + xq) — —2adr, 



xj + x 
and integrating reveals that the asymptotic solution for the deviation 

Ax i (t) = xi (r) + x 

of the instanton from its "vacuum" value —xq decreases exponentially for large To: 

Ax z (r + r) ^ Az 7 (r ) e~ 2Q:T . 



(2.53) 

(2.54) 

(2.55) 
(2.56) 



Moreover, the characteristic time scale 



1 

2^" 



(2.57) 



of the decay of Axj becomes arbitrarily small for large a. In the mechanical analog system, this means that 
to a particle starting from xq nothing much happens for a long time since its velocity remains almost zero. 
However, when it finally approaches the minimum of — V at X = it takes up speed fast, rushes through the 
minimum, and decelerates equally fast near to — xq, spending all the remaining time to creep up fully and 
reach —xq at T — > oo. The "abruptness" of the transition increases with a, the coupling parameter which 
cont rolls the height of the potential barrier. 



3 The instanton should therefore not be confused with a soliton-like particle. There is an important general equivalence, 
however, between instantons and static solitons in the corresponding (field) theory in one additional space dimension. Our 
instanton, for example, is the classical soliton solution (the so-called "kink") of the 1+1 dim. A</> 4 field theory. 
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The (Euclidean) action of the (anti-) instanton solution, which governs the exp onen tial suppression of 
the tunneling amplitude, is easily obtained with the help of the "first integral" Eq. ( 2.4S| ), 



I (YTfl 

Si = Se [xi] —ml drij = m I dxiii = 

-T/2 Jxo x 



i—x 4 

/ dxj (x - Xj) = -amx , 
Jxo •J 



and equal to the action of the anti-instanton. Another way of writing the action, 



(2.58) 



Se [xi] = ml dxiij = I dxy / 2mV (x) 

Jxo J-xo 



dxpE (x) 



(2.59) 



with pe (x) = y/2mV (x) — Ee (recall that Ee = for the instanton solution) shows that the e xpon ential 
suppression factor exp (— Si/H) is nothing but the Gamov factor already encountered in Eq. (1.11). As 
expected, the path-integral approach reproduces the analytically continued WKB result for the tunneling 
amplitude. 

Let us add two comments. First, there exist more saddle points, i.e. solutions of Eq. (2.26), than we have 



discussed above. Those additional solutions do not correspond to particle motion which both starts and ends 
at maxima of — V, however. Hence they do not satisfy the tunneling boundary conditions (2.46,2.47) although 
they might involve barrier penetration. Morevoer, they do not contribute to the SCA for the ground state 
since their action become infinite for T — ► oo. Physically, this is obvious since in such "runaway" solutions 
the particle moves infinitely far from the center of the potential (|x| — > oo) and thereby reaches an infinite 
velocity at T — > oo. Finally, we note that for very low tunneling barriers (i.e. for a — > 0) and correspondingly 
high tunneling rates semiclassical methods may fail. In such cases variational approaches can sometimes help 

0. 



2.2.2 Zero mode 



The constant solutions 1) and 2) above share the full symmetry of the Hamiltonian H (i.e. of V). This is 
not the case, however, for the instanton solution since the latter is localized in imaginary time around tq and 
therefore lacks the continuous time translation invariance of the r- independent H. As a consequence, the 
instanton solutions form a continuous and degenerate one-parameter family whose members are characterized 
by their time-center To- In the saddle-point approximation to the path integral we have to sum (i.e. integrate) 
over the contributions from all tq. In the present section we will see how this can be done. 

To this end, let us go back to our SCA "master formula" ( 2.4C| ) which was derived under the assumption 
A„ > 0, i.e. for positive eigenvalues of the fluctuation operator F . It is not difficult to see that this assumption 
is not valid in the case of the instanton. Since the instanton "spontaneously" breaks the time-translation 
symmetry of V there must be a zero mode in the spectrum of F (for T — > oo), i.e. an eigenfunction 



*o(t) = xl-jg-xi 



■ia 



1 



2 cosh (a(r — To)) 



with the eigenvalue Aq = (which we have normalized to one at T — > oo, where (2.48J) applies 



T/2 /-T/2 

dr [xq (t)] = Sy 1 I dr (mxj) = 1). 

-T/2 J -T/2 



(2.60) 



(2.61) 



In order to v erify that ( 2.60 ) is indeed a zero mode it suffices to take a time derivative of the equation of 
motion (2.26): 



m±i - V (xi) = 



md 2 T - V (xi) xi = -Fxi = 



Fx (t) = 0. 



(2.62) 



The physical origin of this zero mode is quite obvious: it corresponds to an infinitesimal shift of the 
instanton center tq, i.e. an infinitesimal time translation of the instanton solution. Since the resulting, 
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shifted instanton is degenerate with the original one, such a fluctuation costs no action and ( 2.35 ) implies 
that the corresponding eigenvalue of F, i.e. Aq, must be zeroQ 



The instanton solution (2.51) is monotonically decreasing, which implies that the zero mode satisfies 
xq ( r ) < (cf- Eq. ( 2.60 )) and therefore has no node. As a consequence, the zero mode is the (unique) 
eigenfunction of F with the lowest eigenvalue. All the remaining modes have A„ > and for them the 
Gaussian integrals in Eq. (2.38) are well defined. This is not the case, however, for the cq integration 



dc. 



oe" 



2K A C 



dc () 



(2.63) 



which is not Gaussian at all! And neither should it be, since To-shifting fluctuations rj (r) ~ xq (t) are not 
damped even when they are large (i.e. for large cq), due to the To-independence of the instanton action. The 
latter renders the integrand co-independent and thus the integral divergent. 

But this is just the kind of divergence which we should expect anyway from integrating over the infinite 
set of saddle points, i.e. over all To, as discussed above. Indeed, it is easy to show that 



dcg oc drQ 
by comparing the deviations dx from a given path x (r) which are caused by 

small time translations to — * To + dro => dx = — — dro = —xidro, 

dr 



small coefficient shifts cq 



Co 



dco 



dx 



dx 
dc a 



dco = xgdeo 



and 



iidco. 



(2.64) 



(2.65) 



(2.66) 



Equating both deviations (and redefining the sign of dec such that the integrals over Co and To have the same 
limits) gives 

dc = J— dr . (2.67) 

V m 

Thus the integration over Cq can be done exactly and (in the limit T — ► oo) just amounts to summing 
over all one- instanton saddle pointsQ In the presence of a zero mode the expression ( 2.40 ) must therefore 
be replaced by 

-1/2 



Zi(-x ,x ) =Afe" 



Se [xi] 



T (det F[xi]') 



(2.68) 



where the prime at the determinant indicates that Ao is excluded from the product of eigenvalues. Of course, 
the factor T in ( 2. 68] ) becomes infinite in the limit T — ► oo which we will take in the end. This infinity is 
a consequence of t he infinite amount of contributing saddle points and will be cancelled by other infinities 
(see e.g. Eq. (2.17)), leaving the observables perfectly finite as it should be. 



2.3 Fluctuation determinant 

Our next task is to make sense of functional determinants as encountered above, and to calculate 



det F [xj] = det 



-m— + V" (as d ) 



(2.69) 



from the leading fluctuations around the instanton (with the zero-mode contributions removed) explicitly. 
Together with the zero-mode part, this will take care of the O (H) contributions to the SCA. Although 

4 It might be useful to anticipate already at this stage that zero modes of fluctuatation operators, especially in quantum field 
theory, are not always due to continuous symmetries of the underlying dynamics. We will encounter an intruiguing example of 
a topologically induced zero mode below when dealing with Yang-Mills instantons. 

5 Since co measures how far the instanton is collectively (i.e. equally at all times) shifted in (imaginary) time, it is an example 
of a "collective coordinate" . 
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the explicit calculation of ( 2.69 ) is a rather technical exercise, we will go through it in considerable detail. 
One of the reasons is to develop intuition for the calculation of determinants in quantum field theory, as 
they generally occur during quantization of extended classical solutions (besides instantons e.g. solitons, 
monopoles, D-branes etc.). Indeed, all typical features of the calculation generalize rather directly to QCD. 
Hence it can serve as a pedagogical substitute for the actual QCD calculation, which is a tour de force |3l| 
beyond the scope of our introduction to instanton physics. 

Despite their uses, the following two subsections lie somewhat outside of our main focus. Readers not 



interested in the details of such calculations may therefore jump ahead to Section 2.4 without compromising 
their understanding of the remainder of these lectures. 

2.3.1 Warm-up: harmonic oscillator 

To prepare for the calculation of the instanton determinant, we first consider a simpler problem. We look 
for a nontrivial potential V which nevertheless leads to a fluctuation operator F whose eigenvalues can be 
obtained easily and analytically. The most obvious choice is to require V" to be nonzero but x- independent, 
i.e. to work with the potential of the harmonic oscillator 

Vho (x) = -mwV. (2.70) 

The Euclidean analog potential is an inverted parabola. Of course, there is no tunneling in this potential. 



Instead of the tunneling boundary conditions (2.27) we therefore choose 

x cl (±T/2) = 0. (2.71) 

The only solution to the equation of motion (2.26) under the conditions (|2.71) is the constant 



x ho (r) = 0, (2.72) 

corresponding to a particle which stays in "metastable" equilibrium at x = forever. Not surprisingly, its 
Euclidean action SE,ho vanishes and hence the exponential suppression factor exp {—Se/K) of the tunneling 
amplitudes is absent, as it should be. 

To O (H) the Euclidean time evolution matrix element of the oscillator becomes 

-1/2 



Z ho (0, 0) = (o | e - H " oT / fi | o) = N (det F [x ho ]) (2.73) 



where we have redefined F by absorbing a constant factor (det m) into A/": 

F[x ho ] = ~-^+Lu 2 . (2.74) 

Since the classical solution ( [2.72] ) is time-translation invariant, there is no zero mode and the prime on the 
determinant can be omitted. 

We now apply the formal definition of the determinant of an operator as the product of its eigenvalues, 

detF[x ho ] = Y[X ho , n , (2.75) 

n 

with the Xho,n determined by 

F [xho] i n (r) = \ho, n x n (r) . (2.76) 

As intended, the eigenfunctions of F [xho], 

/ y]/KT \ / T17TT \ 

i n (r) = a n sin f -— J + b n cos f -— J , (2.77) 
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and the corresponding eigenvalues 

"ho.n 



\~f) 



■ OJ 



n= 1,2,... 



(2.78) 



can be read off immediately. The coefficients a n) b n are fixed by the normalization ( [2.2C ) and boundary 
conditions ( p.21| ). In order for x n to satisfy the latter, we have to require a n = for n odd, and b n — for 
n even. 

We then formally have 

/ °° r , . 9 
Z ho (0,0)=M[ Yl 






-1/2 



(2.79) 



and since the product is infinite, J\f must be infinite, too, in order to prevent Zho from vanishing identically. 
One way of disentangling the two infinities, or in other words to define the determinant, is to relate it to the 
determinant of the fluctuation operator without potential (to = 0): we simply split the above expression as 



Zho (0,0) =Z free (0,0) J| 



1 



uoT 



-1/2 



(2.80) 



where the first factor corresponds to the propagator of free motion 

-1/2 



Zfree (0, 0) = M 



n 



/n7T\ 2 

VY 



2tt 



rah 

2^rT 



(2.81) 



and the forelast equation above is obtained by inserting the momentum eigenstates (x\p) = exp (—ixp/h) /v27r 



Note that ( 2.81 ) amounts to a calculation of Af from the requirement that it renders the free propagator 
finite and equal to the familiar resul t from basic quantum mechanics. Note, furthermore, that J\f is a property 
of the functional integral measure ( 2.22 ) and does not depend on the dynamics specified in Se [x]. We will 
therefore encounter the same factor below when we return to the instanton problem. 



The second, finite factor in ( 2.80 ) is obtained from the standard formula 

fi(i+£) = 8i ° h< * c) 

71=1 V ' 



TIC 



which can be found e.g. in |25J or by Mathematica. Putting everything together, we finally have 
Z h 



Z ho (0, 0) = M (det F [x ho ]) ' ' 2 = ^/^ [sinh (loT)]- 1/2 



(2.82) 



(2.83) 



Actually, this is the full Euclidean propagator of the harmonic oscillator: since the corresponding path 
integral is Gaussian, the SCA becomes exact! Below, we will be interested in the large-T limit 



Z ho (0,0) 



mhuj 
2tt 



exp (ujT) 



-1/2 



mhuj 



jT/2 



from which we recover, via Eq. ( 2.17 ), the familiar ground-state energy 



(2.84) 



(2.85) 



2.3.2 Instanton determinant 

After having calculated two simpler functional determinants (those for the free particle and the harmonic 
oscillator), we are now prepared to attack the original problem. Our task will be to calculate the fluctuation 
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determinant which appears in the O (h) tunneling amplitude 

_ S E [ x 

Zj(-x ,x ) = Z ho (0,Q)e ~ 



S E [xi] j detF[ Xl ]' 



-1/2 



2irhm 



" 2 det F [x h o 



(2.86) 



around an instanton. Above, we have factored out the root of the harmonic-oscillator determinant in ( 2. 68] ) 
and used Eq. ( [2.73 ) for Zho (0, 0) to eliminate M. This "renormalizes" Zj (— x§,x§) according to the procedure 
of the preceding section and leaves us with the above, finite ratio of determinants. Furthermore, we have 
factored out the lowest-mode contribution to the harmonic-oscillator determinant, A/ lo ,o — > w 2 for T — ► oo, 
to balance the removal of the zero mode from the determinant of F [xj]. 

In order to calculate det F [x{\ , we have to deal with the spectrum of the operator 



where 



V" (xi) 



2a 2 m 



F [xj] = -m^-2 + V" (xi) 



(3x 2 j - Xq) = 2a 2 rn [3tanh 2 a (r - r ) - l] = 2a 2 



cosh a (t — tq) 



(2.87) 



(2- 



(cosh x — sinh x = 1) describes the interaction of Gaussian fluctuations with the smooth, step-like potential 
generated by the instanton. The eigenvalue equation of F [xi] becomes, after absorbing the factor (det m) 
as above into the normalization and specializing to tq = 0, 



4a 2 - 



6a 2 



dr 2 cosh z (ar 

subject to the tunneling boundary conditions 

x x (±T/2) = 



ix (t) = Xx x (r) , 



(2.89) 



(2.90) 



for T — > oo (cf. ( ^.2l[ )). As usual, the eigenvalue spectrum will be discrete for the "bound" states with 
A < 4a 2 and (in the limit T — > oo) continuous for the scattering states where A > 4a 2 . 



Eq. (2.89) is of Schrodinger type and can be solved analytically in terms of hypergeometric functions (cf. 
§25, Problem 4 in pl[ , which makes use of the solution of §23, Problem 5 in the same book). To obtain the 
solutions explicitly, we rewrite (2.89) by means of the standard substitution 



£ = tanh (ar) , 
in the form 

where 



d ad, 2 \ d 

d^ = cosh 2 (ar) df = ° ^ ~ ' & 



^ i-a^ 



dr 2 



(i-e 



d£ 



ii-e 



% 



<% 



dt 



xx (0 = 0, 



e 2 = 4 



A 



(2.91) 



(2.92) 



(2.93) 



The solutions of (2.92) in the interval of interest, £ € [—1,1], are associated Legendre functions (see, e.g. 
||, chapter 8.7): 

M0 = ciP!(0 + <*Q!(0- (2 - 94) 

Let us first consider the case e 2 > 0, corresponding to the discrete levels. The first of the boundary 
conditions ( 2.9C| ), i.e. 

ix (£ = -1) = 0, (2.95) 



requires C2 = 0. The second one, 



ix (£ = !) = 0, 



(2.96) 
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can be nontrivially satisfied only for e = 1, 2 (for which the P| (£) reduce to a polynomials). Thus we have 
two discrete levels, with the eigenvalues 



A o = 0, 
Ai = 3a 2 . 



(2.97) 
(2.98) 



As anticipated, the lowest level Ao is the zero mode which we have already found in Section 2.2.2. The 
corresponding (not normalized) eigenfunction P% (£) = 3 (l — £ 2 ) = 3/ cosh 2 (ar) is proportional to ( 2.60J ) 
and, inserted into (2.89), yields indeed the eigenvalue zero. Both discrete eigenvalues are nonnegative, as 
expected in a repulsive potential. Nevertheless, the modified potential corresponding to the eigenvalues e 2 
has two "bound states" , i.e. the discrete levels found above. 

The calculation of the contributions from the continuum states with e 2 < is more involved but pro- 
vides an instructive example of how to obtain a functional determinant in terms of phase shifts. Since the 
associated Legendre functions are not very convenient to deal with for non-integer e, we rewrite x\ (£) 



e/2_ 



(£), which transforms Eq. (2.92) into 



^-Old-O 



e/2 _ 2 



. e/2-1 



i_e 2 ) 6/ +6(1-0 



=/2 



w(0 = o 



or 



(l"e 2 )^-2 (e+ l)^-( e -2)( £ + 3) 



d_ 



In terms of the new variable u with 



ciu; 



., , « = l-2u, 1-£ 2 =4 U (1- W ), ^ 
this becomes the hypergeometric differential equation 



w (£) = 0. 



dw d 2 w A d 2 w 

' ~dj : ~d^ 2 ^~de 



d 2 d 

u (1 - u) — -= + (e + 1) (1 - 2«) — - (e - s) (e + s + 1) 



■u; (u) = 



(2.99) 
(2.100) 

(2.101) 
(2.102) 



with s = 2. Its solutions for positve e 2 agree with those given above. The general solution for e 2 < reads 



w (u) = c 3 2-Fi 



ifc ife ik 

s,s-\ hi, H ,u 

a a a 



C4U 



ik/a T-i 
2-fl 



ik 

-s, s 4- 1, 1 , w 

a 



(2.103) 



(see, e.g. chapter 9.1 of J25|) where we have defined e = ik/a. The opposite s ign, k — > — fc, interchanges 
incoming and outgoing solutions (see below). For our case of s = 2 the solution ( 2.103| ) can be expressed in 
terms of elementary functions, 



w (u) = c 3 (1 - u)~ ife/a [fc 2 - 3i/fca (1 - 2u) - 2a 2 (1-6(1- u))] 



C4U 



-ik/a 



Qua {—ik + 2a (1 — u)) 
(k + ia) (k + 2ia) 



or, transforming back to £ = 1 — 2u, 



1 1 t\ - lk l a 
w (f ) = c 3 ( -li [fc 2 - 3ifca£ + 2a 2 (2 + 30] 



f'4 



l-£ 



-ik/a 



3a(l-0(-zfc + a(l + Q) 
(ft + ia) (ft + 2ia) 



(2.104) 
(2.105) 

(2.106) 
(2.107) 
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and therefore 



MO 



:i-ey k/{2a) w(o 



C3 






-ik/a 



[fc 2 - 3ika£ + 2a 2 (2 + 3£)] 



C4 



fyr^ 



-ik/a 



1 



3a(l-Q(-tfc + a(l + 0) 
(fc + ia) (fc + 2ia) 



V2vr+e, 

To finally restore the original r- dependence, resubstitute £ = tanh (ar) and use 

-ik/a. 



(VT±L 



,-<*/aln(^) 



together with 



±ar 



In f £3= ) = In ( 1 J C °t { ? T \ ~ Sln !; J '? I = In ( « 
\ 2^/T^l) \ 2 V cosh (ar) =F sinh (ar) j \ 2 



2 y cosh (ar) =p sinh (ar) 

As a result, we can write the full set of solutions with e 2 < as 

i x ( r ) = c 3 e~ lfcT [fc 2 - 3ifcatanh (ar) + 2a 2 (2 + 3 tanh (ar))] 

3a (1 — tanh (ar)) (— iA; + a (1 + tanh (ar))) 



= ±ar- In 2. 



(2.108) 
(2.109) 

(2.110) 
(2.111) 
(2.112) 



c 4 e ifeT 



1 



(k + ia) (k + 2ia) 



(2.113) 



(The constants £3 and £4 are left arbitrary since we do not need to fix the normalization of the i\ (r) or to 
impose the boundary conditions ( 2.90Q for our purposes below.) 

In order to calculate the continuum- mode con tribut ions to the determinant, we have to find and multiply 
the eigenvalues corresponding to the solutions ( 2.113| ). Fortunately, there exists an elegant and efficient 
technique |ll| for doing this , whi ch can be applied to functional determinants in field theory as well. It takes 
advantage of the fact that ( 2.89 ) is a local equation. Thus we can obtain the eigenvalues in the asymptotic 
region \t\ — > 00, where th e pot ential induced by the instanton field vanishes (a benefit of expanding around 
a localized solution) and ( 2.89 ) simplifies to 



dr 2 



Z\ (t) = 0. 



(2.114) 



Here the momentum fc, introduced above, replaces A as the label of the continuum modes. Both are related 
by the "dispersion relation" 

fc 2 = A - 4a 2 



(>o). 



(2.115) 



The solutions of Eq. (2.114) are "plane waves". Since their normalization is fixed (elastic scattering), the 
only effect of the instanton-induced potential can be a fc-dependent phase shift. Moreover, the explicit, 
asymptotic solution reveals that no reflection occurs in this potential. Thus we have 



~ / \ ikr+idu 

x\ (T) oc e 

i\ (t) oc e lkT 



for r 



for 



-00, 



+00. 



(2.116) 
(2.117) 



(Due to the absence of reflection, the linearly independent solutions oc exp (— ikr) do not mix in.) The 
above discussion implies that all the required dynamical information about the eigenvalues is contained in 
the phase shifts of the solutions (2.113) in the limit r — > —00, where tanh (ar) — * —1 and thus 



x\{t) 



c 3 e 



-ikr 



3ika — 2a 



~ ikr 



(1 + ik/a) (2 + ik/a) 
(1- ik/a) (2 -ik/a) 



(2.118) 
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By comparing ( 2.116 ) and ( 2.118 ) we read off the phase shifts 

Yl + ik/a\ (2 + ik/a 
\l-ik/a) \2-ik/a 



Si. = —i In 



(2.119) 



How are these phase shifts related to the continuum spectrum the eigenvalues k, and thus to the deter- 
minant to be calculated? Since the phase s hifts encode the behavior of the scattering solutions at large |r|, 
one would expect the boundary conditions ( 2.90 ) to play a role. And indeed, they are all what is needed to 
establish the relation between the eigenvalues and 5k- To show this, we start from the general solution 



igen,X (t) = Ai\ (t) + Bi\ (-t) 



and impose the boundary conditions (2.9C) to obtain 



Ax\ 



Bix - 



T 



= Ax\[- 



T 



Bix - 



T 



= 0. 



This implies 



and has the solutions 



Xx V 2) _ -ikT-i5 k _ i 1 
«A (?) 



r 



(2.120) 



(2.121) 



(2.122) 



(2.123) 



in terms of the phase shifts 8k ■ Due to the boundary conditions the k n are discrete for finite T and become 
continuous in the limit T — > 00 to be taken at the end. 



It remains to calculate 

det F[xi}' 



Ai n„ = i («:„ 



4a 2 



u- 2 det F [ Xho ] A fto , 2 nn=3 (^ k 2 hon + ^ 



3 n n= i {k 2 n + 4a 
4 



n 



n— 3 



(^ho,i 



4a 2 



(2.124) 



where we have factored out the contribution of the second levels (for T — + 00) 



A, 



A 



fto,2 



3a 2 



^' z 



3 

4 



(2.125) 



to the determinant ratio, using kho,n = nn/T (cf. Eq. ( |2.78 )), and specialized to w = 2a. Since the relevant 
range of k- values in the eigenvalue products of (2.124) will turn out not to contain small k (see below), the 
contribution of the two lowest harmonic-oscillator eigenvalues can be multiplied to the above denominator 
with negligible effect (except for simplifying the ensuing expressions). We are thus led to consider the ratio 



n„ = i (fc 2 



4a 2 



n»=i (Ho, 



4a 2 



= exp \, hi 



n=l 



4a 2 



4a 2 



With 



K = { k ho , n I 



25 k k 



I; » It ti 



T 



(since 8k/T <C 1 in the n-region which contributes non- negligibly for T — ► 00) we have 



In 



4a 2 



4a 2 



In 



1 26kkho,n 



Tk 2 



4a 2 



1 2dkkho,n 



Tki 



4a 



2 ' 



In the continuum limit 



x 1 °° rp nOQ 

^./•(M = — ]>>£;/ (fc„)^- / dkf(k) 

^ — 1 ^ — 1 •'O 



(2.126) 



(2.127) 



(2.128) 



(2.129) 
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\&K — Kho,n+l kho,n 



7r/T) we then obtain 

n„ =1 (kl + 4a 2 ) 
4a 2 



nn =1 (fcL, 



= cxp 



1 f°° dk 2S k k 

■k J k 2 + 4a 2 



(2.130) 



To evaluate the integral in the exponent, we note that 



dk 



In 



4a 2 



2k 



k 2 + Aa 2 



allows us to rewrite 



2kS k 



(l fc 2 +4a 2 



dk—— In 
dk 



k 2 ^ 
4a 2 



r^in[i 

Jo an 



(2.131) 



(2.132) 



(the surface term vanishes since ^ =00 = 0) in terms of the dimensionless variable k = kj (2a). The benefit 
of the partial integration is that the /{-derivative removes the logarithm in the expression ( 2.119J ) for the 
phase shifts from the integrand, 



d8 K __ _.d_, 
dn dn 



1 + 2in 1 + in 



1 — 2in 1 — in 

so that the remaining integral can be done analytically: 

2kS k f 00 . ( 2 



l + K 2 



l + An 2 '- 



,, /c 2 +4a 2 



dn 



Thus 



l + K 2 1 + 4k 2 
det F[xi}' 1 



In [1 + k 2 ] =7rln9. 



- 2 deti^[^ ] 12 



(2.133) 



(2.134) 



(2.135) 



and together with (2.86) and (2.84) we can now assemble our final result for Zj (— xq, xq) at large T 

Zi (-a; ,a;o) = 



xhijj 



jT ' 2 uT 



6S E [ XI ] _M 
e « 



nhm 



(2.136) 



This is the quantum mechanical propagator of the double-well tunneling problem, to O (h) in the semi- 
classical approximation around a single instanton. The exact analogs of both the exponential Gamov factor 
and the preexponential factor y/Sj from the zero mode appear in the one-instanton sector of QCD. 

2.4 Dilute instanton gas 

Up to now, we have concentrated on the saddle points which correspond to singlc-instanton solutions. This 
is not the whole story, however: there are additional (approximate) saddle points which also contribute to 
the semiclassical tunneling amplitude for large T. We will now analyze those "multi-instanton solutions", 
first in the familiar double-well potential and subsequently in its periodic extension (which represents the 
closest quantum-mechanical analog of the semiclassical Yang-Mills vacuum) . 



2.4.1 Double- well potential 

Since the instanton deviates only in a small time interval At — 1/ (2a) appreciably from xo or — xq (cf. 
Eq. (2.57)), and since the overlap between neighboring instantons and anti-instantons is exponentially small 
(cf. (p. 56)), mufti- (anti-)-instanton solutions of (2.26) can be approximately written as a chain (i.e. an 
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ordered superposition) of N alternating instantons and antiinstantons, sufficiently far separated in time by 
the (average) intervall 

A T = § » i-. (2.137) 

Those chains correspond to N tunneling processes, back and forth between both minima of the potential. 
It is intuitively clear that their importance, associated with the frequency of their occurence, increases with 
growing T . 

The approximate A-instanton solutions, composed of single, alternating instantons and anti-instantons 
at times To,fc, can thus be written asn 

N 

x n (t) ='^2x I! j(t -To,k) , (2.139) 



where N must be odd in order to satisfy the boundary conditions ( 2.46 ) and ( 2.47J ). They become exact 



solutions for infinite separations |to t+i — To,fc| — > oo. The (anti-) instanton centers are ordered in Euclidean 

time as 

T T 

— 2 < 7-0,1 <r 0)2 ... < r 0jA r< -. (2.140) 

All multi-instanton solutions have to be included as additional saddle-points in the SCA. We are now going 



to derive the corresponding expression for Z (—xq,Xq) by using the approximate solutions (2.139) instead. 
This simplification goes under the name of "dilute instanton gas approximation" (DIGA). 

The contribution of the approximate A-instanton solution to the path integral is 

Z N ~ N J ' D[i\e- s ^ XN+ ^l h . (2.141) 

We now write the general fluctuation rj (t) as a sum of independent, localized fluctuations rj k (r) around the 
single (anti-) instantons and rjo (r) around the approximately constant pieces x (r) = ±Xo beween them: 

N 

r](T)=Vo(T)+Y,Vk(T). (2.142) 

This implies that 770 (r) can be finite all over [— T/2, T/2] (except at the boundaries) while the rjk (r) are 
time-localized around the fc-th (anti-) instanton. 

Hence the action approximately decomposes into the sum of actions for single (well-separated) (anti-) 
instantons and a "vacuum" piece, 

JV 

Se [xn + rj\ ~ S E [xq + rjo] +^S E [xi + rj k ] (2.143) 

fe=i 

(recall that S [x c i = ±a;o] = and S [xq + rf\ = S [—x$ + 77], due to the symmetry of the potential). This 
formula expresses the (approximate) fact that the (anti-) instantons have too little overlap to interact. 
Furthermore, the path integral measure factorizes into integrals over the localized fluctuations around the 



Eq. (2.139) is to be understood as a shortcut for the more precise (and tedious) expression 



(iV+l)/2 (JV-1J/2 

X N (t) = ^2 X I ( T ~ T 0,2k-l) + ^2 x l( T - T 0,2k) (2.138) 

fc=l fc=l 

which takes explicit care of the facts that the first and last transitions in the chain must be instantons, and that those inbetween 
consist of an anti-instanton and (JV — 3) /2 pairs of an anti-instanton followed by an instanton. 
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instantons and those inbetween[]. As a consequence, Zjy factorizes as 

Z N (-xo, x )~N ( D[r) ]e- SE[xo+vo]/h X f[ M I D[q k ]e' SE[xi+r)k]/h (2.144) 



fc=i 

= Z (x ,x ) [Zi(-Xo,xq)] . (2.145) 

We have written the above product in a sloppy manner. Strictly speaking, the N factors Zi (— xo,xo) 
have different boundary conditions since the endpoints ±xo are (almost) reached at different times in each 
factor, associated with the time-ordering of the (anti-) instanton from which they arise. However, due to 
time translation invariance this makes no difference for the value of the Zi (—xq,Xo) except through the 
zero-mode contribution. Indeed, for the latter we have to integrate over the temporal position r of the 
(anti-) instanton (or, equivalently, over cq), 



S/ ' dr = Zj I dr , (2.146) 



2nhm 



where now the range of the center To,fc of the fc-th (anti-) instanton is restricted by the requirement that it 
occurs after the (k — l)-th, i.e. To,fc-i < To.fc < T/2. The integration over the tq^ thus takes the form 

T/2 j-T/2 /.T/2 rpN 

dro,i / dro.2... / dr , N = -ttt (2.147) 

T/2 Jt 0> i Jr ,N-i iV! 

which results in 

N 

ZiT] 



Z N ~ Z ^ N { . (2.148) 



The above expression determines the ./V-instanton contribution explicitly since we had already calculated Zq 
in (2.84) (with lu = 2a) and since Zi can be immediately obtained from (|2.136) and (2.58): 



Z (±x ,±x )=Ar(det[-d 2 T +u> 2 ])- 1/2 -> (H^.) 1 \-^l\ (2.149) 



~ Zl = 2 J Q ^M e -^P^ = 4 J2a^g e _ |Qmx ■/» (2150) 

V iram V nn 

In order to collect the multi-instanton contributions to Z (xq, —Xq), we have to sum over all odd N (recall 
that odd TV's are required by the boundary conditions (2.46|), (2.47)) and obtain 



( ZjT ) Z f ~ i / \ 

Zdiga(x q ,-x ) = Z J2 N / =Y{ eZlT - e ~ ZlT }= Z v sinh {Zi T )- ( 2 - 151 ) 



N odd 



(Note that the second exponential on the RHS above removes the contributions even in ZjT and doubles 
those odd in ZiT.) An analogous expression for Z^iga (—xq, —xq) is obtained when only contributions from 
even numbers of instantons are summed. Both results can be combined into the expression 



->DIGA 



1 /fi^ 1/2 



(±X , -X ) = t [^L\ e -T/2 j^T T e -i,T} (2152) 

= I f^\ V2 | e -(-/a-^)T T e -(u,/*+*)T j (2.153) 



7 Due to the limited support of the localized fluctuations, they do no longer form a complete set on the whole time interval. 
Any fluctuation x (r) for r € [— S-, -J-] must therefore be described piecewise in bases around every instanton. Since the path 
integral is the product of integrals over the coefficients of those base functions, it factorizes naturally. 
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from which the two lowest energy levels of the system (i.e. those of the ground state and the first excited 
state) can be found as 



En 



Ei 



-H lim — \nZ DIGA (±x ,-x ) = -h lim — 

T — >oc 1 T — >oo 1 



-h lim 



1 r fu 



Zt)T 



htu 



HZt. 



Zr)T 



tko 



hZ, 



(2.154) 

(2.155) 



As expected, the effect of tunneling is to split the degenerate (would-be) ground state energies fno/2 of 
the wave functions |±aro) centered in each of the two minima of the potent ial^ T he corresponding energy 
eigenstates are similarly obtained from the prefactors of the exponential in ( 2.153J ) (cf. Eq. (2.16)). For the 
ground (first excited) state we find the symmetric (antisymmetric) linear combination 



|0) = — {\x ) + \-xv)} 



(2.158) 



\1) = —{\ Xo )-\- Xq )}. 
These are the standard WKB results for tunneling amplitudes, with the typical splitting 



AB~e" 



Sb[» 



(2.159) 



(2.160) 



of the energy levels of the states connected by tunneling, although we have obtained them in the somewhat 
less familiar framework of the imaginary-time path integral. It is remarkable that the SCA works well even 
for the lowest-lying states of the system, i.e. for those with the largest de-Broglie wavelengths A, although 
this A is not small compared to the size of the potential. Note furthermore that, in contrast to the case 
without tunneling and the ground states |±Xo), the above states are eigenstates of parity (under which 
xq <-> — xo). The new ground state ( 2.158 ), in particular, is parity invariant: the artificially broken parity in 
the absence of tunneling is restored. 

Before le aving t his section, we should address a potential concern in summing over the dilute instanton 
gas in Eq. (2.151). Indeed, for an increasing number N of (anti-)instantons in the constant interval T 
the diluteness condition ( 2.137 ) is less and less satisfied, implying that from some large N onwards the 
corresponding terms in the sum will violate this fundamental DIGA requirement. However, Stirling's formula 



/27m for large n shows that the sum (2.151) is dominated by terms with 



ZiT 



0(1) 



and that contributions from larger N are rapidly suppressed. As a consequence, only N with 

N 



T 



<Ce 



(2.161) 



(2.162) 



i.e. with an exponentially small instanton density in the scmiclassical limit (which can be improved at will 
by increasing a since we had found in Eq. ([2.58) that Se [xi] = (4/3) amxg), contribute significantly to the 



8 We note in passing that the energy splitting between the classical would-be ground states could have been obtained from 
the one-instanton approximation alone, i.e. without employing the DIGA. Indeeed, restricting to small T and expanding 



X f 6 Xi 



l--{ Xf \H\x z )+0{T 2 ) 



the propagator is governed by the (time independent) matrix element 

(x f \H\xi) 
which can be calculated in the one-instanton approximation. 



(2.156) 



(2.157) 
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Figure 2.5: A multi-instanton solution in the periodical potential. The abscissa denotes the Euclidean time 
t while the ordinate gives the position variable x. The integer values of x correspond to the minima Xo, n of 
the periodical potential. 



s um (p. 151 ). This ensures that Zdiga is dominated by terms in which the DIGA dilutcncss requirement 
( [2.137D is well satisfied^. 

Our main motivation for the use of instanton methods in the solution of the above tunneling problem was 
that this procedure can be straightforwardly generalized to gauge theories. In order to develop the analogy 
with QCD as far as possible, however, we should still go a step farther in quantum mechanics and study a 
periodical potential with degenerate minima. This will be the subject of the following section. 



2.4.2 Periodic potential 

Let us now consider the periodic extension of the double well potential (which thereby becomes bounded 
from above) . This type of potential most closely resembles the situation which we will encounter below in 
the QCD vacuum. 

In a periodic potential with degenerate minima, instantons and anti-instantons can arbitrarily follow each 
other, starting out at the minimum where the previous one had ended, i.e. connecting adjacent minima i£o,n 
and xo tn ±i (see Fig. 2.5). Generalizing our earlier conventions, we will refer to an instanton (anti-instanton) 
as the segment of the solution which interpolates between neighboring minima to the left (right), i.e. which 
decreases (increases) n. Getting from the minimum with index ni to the one with index n/, corresponding 
to the boundary conditions 

T\ (T\ 

XQ,n t , X[ — \=X0, nf; (2.163) 



x -- 



therefore requires the number of anti-instantons minus the number of instantons, Nj— Nj, to equal n/ — ni 
As a consequence, the semiclassical propagator becomes 



oo oo 



^per yEnfjX ni j — Zjq J > 



ZiT 



N+Nj 



N I= N I= 



N\N r \ 



JNj— Ni— {nf— «j)' 



(2.164) 



With the representation 



Sab = 



d@ i8(a-b) 

„ 2n 



(2.165) 



9 Unfortunately, it will turn out below that in QCD, where we do not have an adjustable parameter analogous to a, the 
situation is much more complex and the corresponding DIGA does not provide a self-consistent approach to the vacuum wave 
functional. 
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of the Kronecker symbol, Z per can be rewritten as 



7 ( \ ~ 7 f^ M ~ 

JO ZlT 



i0(n f -n z ) y^ 
N,=0 



/ - \ N ' / ~ \ N I 

(ZiTe- ie ) °° \ZiTe lS \ 



N T \ 



E 

AT 7 -=0 



JV>! 



(2.166) 



2tt 



dO 



- Z I rL e -^(^f~^) e ZlTe-"'+Z I Te^ 

2tt 



and by using again the expression ([2.149) for Zq we arrive at 



J per 



y^rif ; ■Erii ) — 



ihw 



1/2 



cW 



T/2 i w_ -j$(n f - ni ) ZZiTi 

2tt 



(2.167) 



(2.168) 



mhuj\ 1/2 f 2 * d9 



2n 



-i$(n, f -m) -(oj/2-aZ! cos6)T 



(2.169) 



As before, we can now obtain the low-lying energy levels from the T — ► oo limit. The above expression 
shows that we get, in fact, a continuous "band" of energies parametrized by 9 (since (2.16E) is the sum of 
contributions from the lowest energy levels in (2.14)) 



E (9)=-h lim - 



- 2Zi cos 9 T 



Hlo 



2hZi cos 9. 



(2.170) 

(2.171) 



Of course, it is well known that the energies in a periodic potential form continuous "bands" . Our result 
is just the analog of, e.g., the lowest-lying band of electron states in the periodic potential of a metal. Not 
surprisingly, then, the corresponding eigenstates are the "Floquet-Bloch waves" 



1 



\/27r \ir J 



hjj\ y—^ 



e in9 In) 



(2.172) 



where \n) is the state localized at the n-th minimum of the periodic potential. 

This concludes our discussion of instantons in quantum mechanics. Much more could be said about them 
p6| , their cousins (sometimes called "bounces" ) which mediate tuneneling between nondegenerate minima, 
their connections to large-order perturbation theory in the real-time theory, etc. However, we refrain from 
doing so since we have reached our main objective: to obtain the semiclassical expansion and the ground 
state properties in a potential which mimicks the situation in the QCD vacuum. We will now move on to 
discuss "the real thing" , namely instantons in QCD itself. 
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Chapter 3 



Instantons in QCD 



What is the relevance of the above tunneling discussion for QCD? As anticipated, the answer is that several 
pertinent aspects of the SCAF] generalize to four-dimensional Euclidean (i.e. imaginary-time) Yang-Mills 
theory once we have identified the saddle points of the corresponding functional integral in imaginary time. 
Our first task will therefore be to find the minima of the classical Yang-Mills action. Remarkably, they turn 
out to be determined by the topology of the gauge group and to form the minima of a periodic "potential" 
analogous to the one encountered above. 

The corresponding tunneling solutions, the Yang-Mills instantons, are classical gauge fields with intrigu- 
ing topological properties. The latter induce new physical phenomena without analogy in the quantum 
mechanical examples. Some of these phenomena, including several ones related to the light-quark sector of 
QCD, will be discussed in subsequent sections of this chapter. 

3.1 Vacuum topology and Yang-Mills instantons 

3.1.1 Topology of the Yang-Mills vacuum 

In order to develop some intuition for semiclassical ground-state properties of QCD, let us start as in the 
preceding quantum mechanical examples by searching for the minima of the Euclidean Yang-Mills action. 
Restricting for the moment to the gluon sector, the latter reads 

S[G\ = \J<fix[G%,G$ v \ (3.1) 

= ±Jd i x[E?E? + B?B?}>0 (3.2) 

where the gluon field strength tensor is 

Gfj, v (x) = o^Gv — o v G^ + ig [G^, G v \ = G t , t = — , (3-3) 



2 



and the chromoelectric and -magnetic fields are defined as 



E i ~ C° 4 , Bf — —-£ij k Gj k . (3.4) 



1 We note in passing that there exists a somewhat complementary way of looking at the SCA in quantum field theory, namely 
as an expansion in the number of Feynman-graph loops (which reflects the non-perturbative nature of the SCA from a different 
angle). 
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The corresponding quantum field theory (including the quarks) determines the structure of the QCD 
vacuum, i.e. the unique state of lowest energy on which the Fock space is built. Since QCD is strongly coupled 
and therefore non-perturbative at low energies, we expect the vacuum to be populated by strong fields. (This 
is in contrast to standard QED, for example, where the vacuum contains mostly zero-point fluctuations, i.e. 
weakly interacting electron-positron pairs and photons which can be handled perturbativcly.) 

A measure for the strength of the QCD vacuum fields can be obtained from the trace anomaly, which 
relates the energy density e vac of the vacuum to the phenomenologically known vacuum expectation value 
of the square of the gluon field strength tensor, the so-called "gluon condensate" (renormalized at about 1 
GeV): 

e vac ~ "^(0 |ff 2 G 2 | 0) ~ -~|J « e ^" ( 3 - 5 ) 

This relation shows that the nonperturbative vacuum fields are indeed exceptionally strong: they reduce the 
vacuum energy in a tiny cube of size 10 _15 m by about half a proton mass! As we will see in the remainder 
of this section, part of this reduction is due to tunneling processes mediated by instantons. 

Strong fields can contain a very large number of quanta, and those quanta can become coherent and 
render the corresponding action large compared to H. In other words, such fields behave (semi-) classically 
since quantum fluctuations are of O (H) and thus contribute only relatively small corrections. The above 
reasoning suggests that insight into the vacuum fields of QCD may be gained from a semiclassical perspective. 
In the following, we will explore this perspective while paying special attention to robust and generic features 
which are likely to survive even stronger quantum fluctuations. The most important such features will turn 
out to be "global" , i.e. topological properties of the vacuum fields which are invariant under continuous 
deformations (and therefore in particular under time evolution). 

The first step towards a semiclassical approach to the QCD vacuum (starting at the lowest order of h) 
is to find those classical fields which minimize the (static) Yang-Mills energy or, equivalently, the Euclidean 
action (JO). Such fields are sometimes called "classical vacua". Since the Euclidean action ( |3.1[ ) is non- 
negative (in contrast to its counterpart in Minkowski space) its absolute minima will be the gluon fields with 
zero action. These fields are the "pure gauges" 

G (P9) = —UQTJ\ (3. 6 ) 

9 

where U (x) E SU (3) is an element of the gauge group. It is easy to check explicity that the field strength 
of pure gauges vanishes, 

G% 9) = 0. (3.7) 

At first sight, one might think that these fields would be natural candidates on which to build a semiclassical 
expansion. A moment's reflection shows, however, that none of them could generate an acceptable vacuum 
since they are neither unique nor gauge invariant. 

In fact, it can be easily seen that the G^H fall into an enumerable infinity of topological equivalence 
classes. In order to demonstrate this, we start by choosing the temporal gauge 

Go = (3.8) 

so that the residual gauge transformations (which conserve this gauge condition) are time-independent. 
Moreover, for our purposes it is useful to restrict the residual gauge transformations further by letting them 
approach a constant (chosen to be unity) at spacial infinity, i.e. 

U (x) -> 1 for \x | -> co. (3.9) 

This restriction guarantees that the gauge fields satisfiy definite boundary conditions at the surface of a large 
box (which should not affect the local physics inside). 



The effect of the condition (3.9) is that, from the point of view of gauge transformations, the space R 



3 
space 



is effectively compactified to Sj e (i.e. to a 3-dimensional sphere of infinite radius), i.e. different boundary 
points at |x| — > oo cannot be distinguished by U (x) and thus can be identified. This is completely analogous 
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Figure 3.1: Schematic view of the generalized stereographic projection which identifies compactified R 3 with 
S 3 . (Figure from hep-th/0010225 , courtesy of Falk Bruckmann.) 



to the compactification R 2 — > S 2 of the complex plane by Riemann's stereographic projection, as shown in 



Fig. 3.1. As a consequence, any residual gauge transformation U (x) defines a map 



S: 



space 



ST/ (3) 



color ' 



(3.10) 



For the purpose of the following topological (homotopy) classification, this map can be furthermore restricted 
to a SU (2) subgroup of the gauge group, since according to a powerful theorem by Raoul Bott [^7| only this 
subgroup will be "topologically active" in our discussion. Moreover, the quarternion representation of any 

UgSU(2), 



U = u + iu a r a 



with 



real and 



u a u a = 1, 



(3.11) 



(in terms of the Pauli matrices r a ) sh ows that SU (2) can be mapped onto the 3-sphere S^ roup (which is its 
group manifold). Thus the map (3.1C) is topologically (in the sense of homotopy theory) equivalent to 



space group 



(3.12) 



i.e. any residual gauge transformation maps the spacial 3-sphere into the group 3-sphere. In general, two 
such mappings U (x) cannot be continuously deformed into each other, which implies that they fall into 
different homotopy classes. The same holds for spheres of any dimension n, i.e. for any map S n — > S n . 

To get a visual understanding of this topological classification, it is thus sufficient to consider just the 
simplest case 

(3.13) 



9 1 

w sp 



9 T 



of maps between two circles. Such mappings can equivalently be though t of as phase fields U (x) = e*^(*' <E 
U (1) defined on a circle with coordinate angle \ G [0, 2tv]. In Fig. 3.2 various maps of this type are 
graphically represented by drawing the domain space circle (which can be imagined as a rubber band) in 
contact with the target space (the dotted circle which represents the phase of the field) such that the mapping 
occurs between those points on the circles which touch each other. Note that we consider only continuous 
maps which implies that every point of the target circle must somewhere touch the domain circle. The 
examples drawn in Fig. ft. 2| illustrate the topological properties of such maps. We start with Fig. |3.2| a where 
the whole domain space shrinks to a point (which we draw as a tiny circle with radius r — > implied) and 
touches the target space at one and the same point. Clearly this graph represents the constant map: all 
points of "space" are mapped into one and the same point in the field or target space. 



Let us now consider Fig. 3.2b where the mapping (field) becomes space-dependent. We define that fields 
which can be continuously deformed into each other (i.e. such that the all points of the domain space stay in 
contact with points on the domain circle) belong to the same "homotopy class". Now the map of Fig. 3.2b 
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a) n = 



b) n = 




c) n = +1 



d) n = -1 



e) n = 2 



Figure 3.2: Some characteristic examples of the mapping S 1 
dotted circle is the range of values of 0, while the outer, solid ' 



'circles' 



exp (ix) i— ► exp (icf> (%))■ The inner, 
symbolize the domain of %. 
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can clearly be continuously deformed into the identity map and thus belongs to the same homotopy class. 



This is in distinct contrast to the remaining maps of Fig. |3.2| . In Fig. ^c the domain space wraps once 
counterclockwise around the target circle and therefore cannot be continuously shrunk to a point. Therefore, 
the corresponding map lies in a d isjoint homotopy class with which we will associate a "winding number" 
n = +1. The mapping of Fig. p^ d wraps clockwisqj around the target circle and belongs to the homotopy 
class with winding number n = — 1. (The trivial class containing the identity map obviously has n = 0.) 



The next graph, Fig. 3.22, shows a map with winding number n = +2 where the target space wraps 
twice counterclockwise around the domain space. Now it should be clear how this classification proceeds to 
higher winding numbers. We have thus made plausible that the maps ( 3.13|) fall into an enumerable infinity 



of disjoint homotopy classes characterized by an integer winding number n £ Z . If one additionally defines 
a composition law for such maps by concatenation (under which the winding numbers of the maps to be 
composed simply add) the homotopy classes become elements of a group, the so-called homotopy group 

7Ti (S l ) = Z. (3.14) 

Here the subscript of 7r indicates the dimension of the domain sphere (in the present case equal to one), 
and its argument denotes the target space (here S 1 ). Now, according to what we have said before this 
result generalizes to the homotopy groups for mappings between spheres of dimension d, tt^ (S d ) = Z, which 



includes the maps (3.12) which interest us in the context of the QCD vacuum. We thus conclude that 

tt 3 (S 3 ) = Z. (3.15) 



Let us summarize what we have learned so far: the pure-gauge fields ( p.6[ ), constructed from gauge 
transformations U^ n ' (x) which satisfy the boundary conditions ( |3.9| ), minimize the Euclidean QCD action 
and fall into disjoint homotopy classes characterized by an integer winding number n which derives from the 
topological properties of the CA n ) (x). Thus the topology of the gauge group on the compactified space S 3 
induces the semiclassical vacuum structure with its periodic "potential" (action). 

Note that the above arguments do not prevent us from transforming a pure-gauge field of class n, i.e. 

Gj,") = — U {n) 3^U {n) \ (3.16) 

by continuous deformation into one of class m ^ n, despite the fact that we have just shown that we cannot 
continuously deform U^ n ' into U^ m K Indeed, the latter implies only that we cannot go from G^ to G^ 



without leaving pure gauge, i.e. we cannot keep the field always in the form (3.6). This also means that 



we cannot stay in the sector of fields with zero action (recall that the pure gauges are the only gluon fields 
which minimize the Euclidean action). In other words, while continuously deforming GjT into G™ with 
n/mwe necessarily encounter field configurations with non-minimal action, Se > 0: each topological class 
corresponds to an absolute action minimum, and those minima are separated by a finite-action barrier (the 



so-called sphaleron barrier). This situation is depicted in Fig. 3.3 



3.1.2 Yang-Mills instanton solution 

The rea der h as probably already noted the analogy between the d egenerate classical ground state found in 
Section |3.1.l| and the periodic potential discussed in Section 2.4.2 . In both cases, the Euclidean action has 



an enumerable infinity of degenerate minima. And as in the quantum mechanical examples, the degeneracy 
between the classical "would-be" or "candidate" vacua of QCD is lifted by tunneling of gluons through the 
finite action barriers. 

Thus we know already, at least in principle, how to obtain the semiclassical approximation to the tun- 
neling amplitude between the classical minima G)% : it is generated by a saddle-point approximation to the 



2 Note that the direction in which a point on the curve proceeds if the curve parameter increases does not change under 
continuous deformations. 
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Figure 3.3: The classical, Euclidean QCD Lagrangian for different classes of gauge field configurations. 
The degenerate absolute minima correspond to pure gauges ^4( n ),/i (— G> ) with winding number n. An 
instanton trajectory, interpolating between Am and -<4(2)> is also shown. 



Euclidean functional integral around stationary (and approximately stationary) field configurations, i.e. the 
solutions of the Euclidean Yang-Mills equation 



0(jTj/ 



(3.17) 



with the boundary conditions 



G»(3,T = -oo) = GV>(3), 

G„(x,T = +oo) = G^(x). 



(3.18) 
(3.19) 



According to the above boundary conditions, the solutions start out at T — ► — oo as pure-gauge fields 
constructed from a gauge transformation with winding number n, and they end up in a pure gauge associated 
with winding number matT-> +oo. Moreover, they do by "paying" the minimal amount of action possible. 
These solutions are the QCD-instantons, and they can be found analytically [[y, Q. For m = n + 1, the 
instanton (in a non-singular gauge) has the explicit form 



Gi I] (*) = - 



^ ^a\iv \&v %v) t 



g (x - z y + p 2 



(3.20) 



An instanton with n = 1, 



2 is drawn in Fig. 3.3 



A couple of remarkable features can be read off directly from the instanton solution (3.2C): first, its "spin" 



(i.e. the Lorentz vector index) is coupled to the color orientation by the 't Hooft symbol rj atll/ . Second, its 
nonpcrturbative character (as expected in the context of tunneling) reveals itself in the diverging weak- 
coupling limit. (The 1/g behavior is common to all solutions of the classical Yang-Mills equation (3.17), as 
can be seen by rescaling the gluon field.) Finally, the solution makes explicit that the tunneling process is 
localized in a region of size p in space and time, around its center z. 

By using the definition 

"Ha^v — SanS^u — S av S4 l _ l + £ at iv (3-21) 



of the 't Hooft symbol, it is easy to verify that the instanton (3.2C) indeed satisfies the boundary conditions 
(3.18), (3.19) with n = 0. Thus, QCD instantons describe localized, flash-like rearrangements of the vacuum 
which mediate tunneling processes between topologically distinct pure-gauge sectors. 
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In analogy with the quantum-mechanical example of the periodic potential, those tunneling processes lift 
the degeneracy between the pure-gauge "would-be" vacua \n). Instead, superpositions of all those states, 
the "theta vacua" 

lO^^e^ln), (3.22) 

n 

become the gauge-invariant eigenstates. Although part of our above discussion (including the tunneling inter- 
pretation) depended on the gauge choice (|3.8[), the ensuing vacuum structure is therefore gauge-independent, 



too. We will further elaborate on the 9 vacuum structure in Section 3.4. Here we only note that the analogy 
with the Bloch waves of the quantum-mechanical periodic potential is limited. While the symmetry gener- 
ators in the periodic potential (i.e. the generators of periodic translations) are physical, they just generate 
gauge transformations of nontrivial topology in the Yang-Mills case. Moreover, the analog of the Bloch 
momentum, the angle 9, classifies different theories and cannot be changed inside a given "6- world". At 
present, the value of 9 can only be determined phcnomenologi cally. This task is made easier by the fact that 
finite values of 9 induce CP-violating amplitudes (see Section 3.4.2J ). 



Evidence from instanton- liquid vacuum models suggests that approximate saddle points, namely super- 
positions of instantons and anti-instantons |38], E|, dominate the SCA to the generating functional of QCD. 
Such superpositions are approximate solutions of the Yang-Mills equation if the typical separation between 
the (anti-) instantons is much larger than their average size, a condition which seems to hold in the QCD 
vacuum (see below). In contrast to the quantum mechanical examples of Sections 2.2 and 2.4, however, the 
classical interactions between Yang-Mills instantons are of dipole type (at large separation) and therefore of 
much longer range than the exponentially suppressed overlaps between quantum mechanical instantons (cf. 



Eq. (2.56)). As a consequence, the dilute instanton gas approximation fails in QCD. The strong correlations 
among QCD-instantons generate an ensemble which probably resembles more a liquid than a gas. 

The size distribution n (p) of instantons in the vacuum is currently under active investigation, mainly on 
the lattice ||. The obtained results are not yet fully consistent with each other but conform more or less 
to the older results of instanton vacuum models H . (Some calculations have led to a larger density which, 
however, is difficult to measure reliably on the lattice.) For the lowest moments of the distribution, the 
standard values are 

dpn{p) ~ 1 fnT 4 , (3.23) 



p=- I dppn{p)~\im. (3.24) 

n J 3 

Thus, in the QCD vacuum tunneling happens on average in about 5% of spacetime, and it happens very 
rapidly: the action barrier is penetrated almost instantaneously (T tunne i ~ 0.3 x 10 _15 m ~ 10 _24 s), therefore 
the name "instanton". 

3.2 Including quarks 

Some of the most striking effects of QCD instantons are associated with the existence of light quarks, and 
not even a terse introduction to instanton physics like the present one would be complete without mentioning 
at least a few prominent examples. In the following two sections we will therefore try to give a flavor of the 
most important quark-related instanton effects. 

3.2.1 The chiral anomaly 

Behind many essential instanton effects in the light-quark sector lurks the axial anomaly, i.e. the fact that 
the (flavor-singlet) axial quark current, which is conserved in the classical theory, ceases to be so at the 
quantum level. In the following we will derive this anomaly and some of its implications^. Also from a more 

3 A transparent introduction to QCD anomalies (and low-energy theorems) can be found in [EgJ. 
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general point of view, the discussion of anomalies fits well into our context since it exhibits a particular 
strength of the semiclassical approximation: all known anomalies are fully determined by one-loop Feynman 
diagrams, and as such appear exclusively at O (h) of the SCA. 

To set the stage for the following discussion, let us switch on the fermionic part of the QCD action, 

S q [q, q,G} = J d 4 xq [t 7/i (0„ + igG^)] q, (3.25) 

which couples the previously considered gluon fields to light quarks (we omit the small quark mass term). The 
Nocther procedure and the Dirac equation imply that this action has a classically conserved axial current, 
i.e. 

dp, (qhp,i5q) = o. (3.26) 

(Including the factor i is a convenient convention in Euclidean spacetime.) Now let us compute the quantum 
expectation value of this current in the background of a fixed, classical gauge fieldn G M , i.e. 

J 5 , M (x\G) = Z- 1 [G] J VqVq {qi W q) e -S,[«,?,Gfl (3 . 27) 

where 

Z[G]= { VqVqe~ s " [q ^ G] . (3.28) 



In order to perform the standard Grassmann integration over the quark fields in (3.27), we introduce 
Grassmann-valued sources 77, 77 and write 



J 5l „ (x\G) = Z- 1 [G] (J-i^l.} Jv q Vqe-I di ^ v «+™+^ 



(3.29) 

77.^—0 



where we have defined the Euclidean Dirac operator 

V = i 7M % + igGp (x)] . (3.30) 

With the help of its inverse, the Dirac propagator S (x, y\G) in the background field G M , defined by 

[*7/i (ft,,/. + *<?G M (a;))] S (x, y\G) = 5 A (x - y) , (3.31) 

we can rewrite the exponent in the integrand of (3.29J) as 



qDq + rjq + qr\ = (q + r]S) V (q + Sr]) — rjSrj (3.32) 

and, after redefining q — > q' = q + 5*77 (this is a constant shift of the functional variable q which leaves the 
measure invariant), we obtain 



h 



VqVqe~ $ d x li' D '!+v<}+<iv] _ 2, \G\ e^ d xd v^WSix^GWy)] _ (3.33) 

Inserting this identity into Eq. (3. 29]), we find 



J 5 , M (x\G) = -N f tr ltC {ii^S (x, x\G)} (3.34) 

(Nf is the number of quark flavors; the trace is over Dirac and color indices). This expression is still formal: 
the singular coincidence limit of the propagator has to be regularized, and we will do so below. In terms of 



Feynman graphs, Eq. ( p.34|) states that the expectation value of the axial current can be calculated from 
a closed quark loop in the background of a (classical) gluon field G, with one insertion of the axial vertex 

7m75- 

4 Note that this is the first step towards a full quantum calculation, where we would afterwards integrate over all classical 
gauge fields. In order to derive the anomaly, however, it is sufficient to quantize just the fermion sector. 
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For the further evaluation of ( 3.34 ) we employ the spectral representation of S (x,y\G) in terms of the 
normalized eigenfunctions %p n of the Dirac operator T>, 



Vlp n (x) = A„Vn (x) , 

(T> is hermitean (7+ — 7^ and G+ = G M ) so that its spectrum is real) which reads 



S(x,y\G) 



E j jx) tpn (y) 



(3.35) 



(3.36) 



(For the moment we do not worry about vanishing eigenvalues, although they will play a prominent role 
later onf].) In order to regularize this infinite sum we introduce a cutoff function 



U (a?o 



-eK 



which suppresses the contributions from large |A„| and goes to unity in the limit e 
to take in the end. Inserted into ( 3.34 ) this yields 



J 5 ,„ (x\G) = -N t £ ^{x)i^n{x) e _ sXl 



(3.37) 
which we are going 

(3.38) 



Let us now obtain the divergence of this current by noting that ( |3.35 ) implies 

dp (lpnh^l5lpn) = 2A n (^„75-0n) ) 

and leads to 



dp J 5, M ( X \ G ) 
where 



-27V/ J2 V-n75e" 



-:K 



i>n 



-2N f J2 i>n75e~ eV '' Vn = 2N f tr^ c (x 



75 e 



-cC 



V 2 = -(dp+igGp) -la^G^ 



(3.39) 



(3.40) 



(3.41) 



(<r»v = (i/2)[7/i,7v])- 

We can now expand the exponential in powers of the classical background field G. In Feynman-diagram 
language, this corresponds to expanding the quark loop in the number of G insertions. Each additional G 
insertion implies an additional quark propagator in the loop and therefore makes the loop integral converge 
faster at high momenta. Thus there can be at most a few UV-divergent terms, associated with the lowest 
orders of G. These terms require regularization, and it is at this point where the classical axial symmetry 
gets broken (if we insist on keeping the vector current conserved) and the anomaly emerges. 

The above reasoning is reflected in our calculation by the fact that only the leading nonvanishing term 
in the G expansion of ( |3.40 ) will survive the limit e — > which we are going to take in the end (since higher- 
order terms in the expansion of the exponential (3.40) also introduce higher orders of e). We therefore 
just have to identify and calculate this single contribution. Since all the terms from the first part of (3.41) 
vanish due to ir 7 (75) = 0, and since those from the first-order contribution of the second term vanish due to 
tr 7 (750^1,) = 0, the first nonvanishing contribution originates from the second-order piece of the c^-term, 
i.e. 



tr 



-f,c 



75 e 



- £ B k 



^d 2 



S 



-tr 1 (75CT ^o pa ) tr c {GpyGpu) + O (e 3 



(3.42) 



which corresponds to a quark loop with two (vector) interactions with the background gluon field (cf. Fig. 
3j). With 



5 Alternatively, one could regulariz e the divergence they cause in (3.36) by introducing a small quark mass. 

6 Note that the G-dependence of (3.40) originates solely from the exponential since the implicit G-dependence of the ^ n 

eT> 2 ' 



cancels in the sum. (This is because the matrix element (x 75c 
may decide to evaluate it.) 



\ x ) is independent of the basis of functions in which one 
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Figure 3.4: The triangle graph correponding to the UV-divergent quark loop which generates the chiral 
anomaly. The bigger blob denotes the insertion of the axial- vector vertex 7^75, the smaller one the insertion 
of the vector vertex 7^ to which the gluons are attached. 



and 



TTry X^I^G iiv@ pa) ^pupaj 

)= fd 4 p(x\p)e ed2 (p\x) = 



*-* [iv ~ £p„vpa ^ pa 



d 4 P 



1 



(2-r 



16tt 2 £ 2 



(3.43) 



(3.44) 



(note that this term contains the UV singularity for e — > 0), where we have used (x\p) = exp (—ip^x^) / (2-7r) , 
we finally obtain (in the limit e — > 0) 



d^J^ (x\G) = ^ftr c (G^G^) 



(3.45) 



This famous result (which can be shown to be independent of the regularization procedure) expresses 
the gist of the "axial anomaly" : although the axial current is conserved in classical chromodynamics, the 
corresponding symmetry is broken at the quantum level, i.e. in QCD, due to the necessity of regularization 
and renormalizationFL 



3.2.2 Quark zero modes and index theorem 

The existence of the axial anomaly has crucial implications for fermions in the background of instanton fields. 
In order to exhibit them, let us take a closer look at the Dirac spectrum 



Vi> n (x) = \ n ip n (x) 
introduced in the last section. Since the (hermitean) Dirac operator anticommutes with 75, 

{©,75} = 0, 

each (normalized) eigenfunction tp n with A rl 7^ implies the existence of another eigenfunction 

ip- n (x) = 7 5 V>« (x) 
with eigenvalue 



A_, 



-A. 



(n > 0). 



(3.46) 

(3.47) 

(3.48) 
(3.49) 



In other words, nonvanishing eigenvalues appear in pairs with opposite signs. The eigenfunctions ipo,k, 
corresponding to the remaining eigenvalues Afc = 0, are called zero modes. The zero modes can be made to 
diagonalize 75 (due to (3.47)) by transforming to the chiral basis 



^0±,fc = r ( 1:t 75)^0,fc 



(3.50) 



7 Readers who are familiar with the original derivation of the anomaly in terms of explicit Feynman graphs will have noted 



that we just calculated the celebrated triangle graph Fig. 3.4 in a somewhat implicit, functional fashion. The functional 
approach allowed us to get the desired result without introducing the machinery of perturbative quantum field theory (i.e. 
Feynman rules). 
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so that 

75^o±,fc = ±Vo±,k- (3-51) 

Since, finally, the QCD Dirac operator is flavor-independent, i.e. 

[D,T a }=0, (3.52) 

(where T a are the generators of the flavor group) each zero mode comes in Nf copies (one of each flavor). 

The zero modes V'o,± & re normalizable and inherit several characteristic properties of the instanton, 
including the localization in space and time and the coupling of spin and color. In the so-called singular 
gauge, their explicit form |3l[] is 

Vq,±{x) = - 7-0- 2W2 - u. 3.53) 

(The spin-color coupling matrix u is defined by (a + f) u = 0, and r = x — Xo.) For our subsequent discussion 
only qualitative properties of the zero modes will matter, though. 

The two main characteristics of the Dirac spectrum established above - paired eigenvalues of opposite 
sign and chiral zero modes - have profound consequences. One of those becomes explicit when one integrates 



the anomaly equation (3.45) over spacetime and uses Eq. (3.40) (multiplied by a factor 1/ (27V/)) to write 

Q = y|^ / d4xtr c (G^G^) (3.54) 

= ~X! e " A " / ^^nlb^n- (3-55) 

Above we have defined the "topological charge" Q whose meaning will be clarified below. Since the eigen- 
functions ipn an d 75^ of V have different eigenvalues (A„ and — X n ) and are thus orthogonal to each other, 
the integral in (|3.55 ) vanishes for all n with A„ 7^ 0, leaving us with an integral over the zero modes only: 






Q = -J~] d 4 xipo, k j 5 ip , k (3.56) 



n - r n + r 

= X/ / ^^O-.mV'O-.m - ^ / d 4 X^o +tm tpO+,rn (3.57) 

where n± is the number of right- (left-) handed zero modes per flavor in the given background gluon field. 
Since the eigenfunctions of T> (including the zero modes) are normalized, we finally obtain 

Q = n--n+. (3.58) 

The above formula, relating a property Q [G] of the gauge field to the number of unpaired quark zero 
modes of the Dirac operator ( |3.35 ) in its background, is a special case of the celebrated Atiyah-Singer index 



theorem ]29|, [3JJ for the Euclidean Dirac operator. It implies that Q must be integer and thus cannot change 
under smooth variations of the background gluon field. Such a behavior is characteristic for topological 
properties of the gluon field and prompts us to have a closer look at those in the next section. 

Before doing so, however, we should mention a few phcnomenological consequences of the instanton- 



induced zero modes in the quark spectrum. The integrated anomaly equation (see Eq. (3.62) below) reveals 
that the Q = 1 instanton is accompanied by the appearance of 2Nf units of axial charge. This can happen 
either by creating Nf left-handed quarks and annihilating Nf right-handed antiquarks in zero-mode states or, 
in the alternative (crossed) time-ordering, by helicity flips of quarks propagating in the zero-mode state. Such 
processes are described by non-local 2A^-quark interactions known as 't Hooft vertices [plj. They generate 
the dominant instanton effects in the light-quark sector and crucial interactions between instantons. As 
a reflection of the anomaly, the 't Hooft vertices manifestly break the axial U (1) symmetry of the QCD 
Lagrangian with m q = and thereby resolve (at least qualitatively) the so-called "E/(l) problem", i.e. they 
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explain why the rf meson has almost twice the mass of the r\ meson and cannot be considered a (quasi) 
Goldstone boson. More generally, the quark interactions mediated by 't Hooft vertices have a characteristic 
channel dependence which makes them, e.g., strongly attractive in spin-0 quark-antiquark channels but mute 
in vector and axial-vector channels. These patterns are reflected in the instanton contributions to hadronic 
amplitudes. We will come back to this issue in the final lecture. 

3.3 Instanton topology 

So far, we have established two important implications of the axial anomaly: first, in order for the anomaly 
to have physical significance there must be gauge fields with non- vanishing topological charge Q, and second, 
such fields imply the existence of unpaired zero modes in the spectrum of the associated Dirac operator. 
Our next tasks will be to better understand the topological charge, the gluon fields which carry it, and the 
impact of both on the vacuum structure. 

3.3.1 Topological charge 



We start by noting that the integrand of the topological charge Q (cf. (3.54)) can be written as a total 
derivative 



9 



16tt 2 



:tr, 



I {- T fAl/(- T lll' J 



d^ 



(3.59) 



(we do not need the explicit expression for the "Chern-Simons current" K^ here), which implies that only 
fields with nontrivial behavior at the spacetime boundary x 2 — > 00 (typical for topologically active fields) can 
carry a finite Q. Furthermore, the index theorem ( 3.58| ) implies that Q can only take integer values and cannot 
therefore change under continuous deformations of the gluon background field. Thus, besides the homotopy 
classification of the pure gauges encountered in Section |3.1.l| , gluon fields carry another topological property 
specified by their Q-values. As we have seen above, fields with nonvanishing Q "activate" the anomaly by 
generating quark zero-modes and are therefore, via Eq. ( |3.45| ), associated with a non-conservation of the 
axial charge q§ (t) = J d xJs^ (x\G): 



Q = j^ / d 4 xtr c [G^G^j = — - / d A xd^ J 5)/ 



(x\G) 



(3.60) 



2N f 



d 6 x / dr 



diJ 5A (x\G)+VJ 5 (x\G) 



2N f 



d 3 xJ 5A (x\G) 



T— — OO 



(3.61) 



2N f 



Ag 5 



(3.62) 



(3.63) 



where 

A<?5 = <75 (t = °°) ~ 15 (t = -00) 
is the amount of axial charge created by the gluon field with topological charge Q. 

We will now show that an integer value of Q can be assigned to any gluon field with finite (Euclidean) 
action (only those are relevant for the SCA), and that this value has a transparent physical interpretation. 
For the action (3.1) to be finite, the field strength G^ v has to be square integrable, i.e. it must vanish 
towards the boundary S 3 (x 2 = 00) of Euclidean spacetime, i.e. for |x| — ► 00, faster than 1/x 2 . Now we know 



from Section 3.1.1 that the Euclidean action has its absolute minimum at G^ v = 0, i.e. for pure gauges. As 



(3.64) 



a consequence, finite-action fields satisfy the boundary condition 

lim G fi = —U^(x)d^ n) - 1 (x) 

x 2 — >oc CJ 

where U^ n ' is an element of the gauge group Q with winding number n. (Not surprisingly, instantons are 
of this type (cf. ( |3.18| ), (3.19)) as they interpolate between pure gauges of different winding number for 

T — ► ±00.) 
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The essential consequence of Eq. ( 3.64 ) is that it implies a time-independent topological classification 
of all finite-action gauge fields. This classification is based on (but not e qual to) the ho moto py analysis of 
the gauge group elements U^ which we went through in Section 3.1.1 . Indeed, Eq. ( 3.64 ) implies that 
any finite-action gauge field defines a map S 3 (x 2 — oo) — > Q, and such maps fall into disjoint homotopy 
classes according to the third homotopy gro upp| 7r 3 (g) of Q. For Q = SU(N) with N > 2, which includes 
QCD, we have ir 3 (SU(N)) — Z as in Section 3.1.1 while, e.g., (non-compact) QED has 7r 3 ([/(l)) = 1. As a 
consequence, all finite-action gluon fields in Euclidean QCD fall into topologically distinct equivalence classes 
labeled by the topological charge Q (which is technically the "Pontryagin index" of the gauge field )|j. 

For instantons, the topological charge Q has a particularly transparent meaning: it is just the difference 
between the winding numbers n and m of the pure gauges between which the instanton interpolates, 



= m — n. 



(3.65) 



Instead of proving this expression in general (which would not be difficult), let us check it for the explicit 
instanton solution (3.20) found above. We first calculate the classical instanton action (from which we know 
already that it is finite and that it is minimal in its Q-sector) by plugging the field strength 



Gi^ v (x) 



-Ap 2 



^ja^iu^a 



9 [(x-xq) + P 2 



(obtained by inserting the expression (3.20) for the instanton into the Yang-Mills field strength tensor (|3 
into the Yang-Mills action (3.1). The result is 



(3.66) 



)) 



St 



d A xtr c {Gi tllv Gi 



///'/ 



8^ 

«2 ' 



(3.67) 



We note that the classical action is additive for multi-instantons and multi-anti- instantons, Sq = \Q\ x Si 
(but not for combinations of both), and that it does not depend on the instanton's collective coordinates 
Zp, U and p. The latter is a direct consequence of the translation, gauge, and scale invariance of classical 
Yang-Mills theory. At the quantum level, the scale invariance gets broken by the trace anomaly (yet another 
variety of anomaly) and the effective quantum action becomes p-dependent. This generates a p-dependcnt 
weight n {p) in functional integrals which can be identified with the instanton size density introduced at the 
end of Section 3.1.2. 



Having calculated the instanton action, we can easily check whether the instanton has indeed Q = 1. 
In fact, this becomes trivial once we have established another essential instanton property, namely the 
self-duality of its field strength. We now turn to this issue. 

3.3.2 Self duality 



Inspection of the instanton's field strength (3.66) reveals that it depends on the Lorentz indices only via the 
't Hooft symbol i] ailv . Moreover, the expression (3.21) for r\ a y,v shows that it is self-dual (recall the definition 
Q), i-e. 

which directly carries over to the field strength of the instanton: 

Gi 4LU = Gi^y. (3.69) 

Note that self-duality is a gauge invariant property. It implies Ef — Bf , i.e. the equality of the chromoelectric 
and chromomagnetic fields (3.4), and yields immediately 



Q' = ^2Jd 4 -tr c {Gi^G i^) = JL; J 



, , d i xtr c (Gi^ u Gi. l _ llJ ) = 

107T" ' 



JL 

8tt : 



■Si = 1. 



(3.70) 



8 In mathematical terms, 713 (G) classifies principal Q fibre bundles over S 4 [B0|. 

9 At this point, it might be usefu l to emphasize the differences between the topological classification discussed here and that 
of the pure gauges in section B.l.lL The winding number n of the latter classifies maps from compactified space Sf pace (on 
which static gauge transformations with unit-value at the boundary are defined) into the gauge group, while Q characterizes 
maps from the boundary S 3 of Euclidean spacetime into the gauge group. 
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The anti-instanton is anti-selfdual, i.e. Gj „ v = —Gj„ v , and thus has Qj = — 1. (Anti-) Self-duality 
is a mathemati cally very powerful property. Self-dual fields automatically satisfy, for insta nce, t he Yang- 
Mills equation ( p,17| ). This is immediately obvious when combining the self-duality equation ( 3.69 ) with the 
Bianchi identity 

£>,A„ = (3.71) 



which holds for any gauge field, due to the form (3.3) of the field strength tensor. In practice, it is often 
advantageous to deal with Eq. (3.69), instead of with the Yang-Mills equation itself, since it is of first 
order. As a matter of fact, the original instanton solution was found in [p| by solving (3.69|). Similarly, in the 



quantum mechanical example of Section 2.2.1 we derived the instanton solution by solving the corresponding 
first-order equation (2.48). Self-duality has even more profound mathematical implications, e.g. in Donaldson 
theory. It can be shown, incidentally, that all minima of the Euclidean Yang-Mills action with Se < oo are 
(anti-) self-dual, and that they correspond to (mufti-) instanton solutions fll3fl . 



The index theorem ( 3.58 ) implies that there must be one unpaired, left-handed zero-mode per flavor in a 



gluon field with Q = 1. This result can be sharpened for (anti-) selfdual fields with G^ = ±Gfi V . In order 
to decouple left- and right-handed modes we consider the eigenvalue equation of the iterated Dirac operator 
(Cf. Eq. (plb), 



V 2 i, n 



(d^ + igG^) - '-a^G^ 



Ipn = A„V« 



(3.72) 



For (anti-) selfdual fields the second term in the square bracket can be rewritten with the help of the identity 



' [If 



-7507W 



as 



OwGffi = Tv^Glfjn = a^Gffi- (1 T 7s) ■ 



Restricted to quark zero modes and to (anti-) selfdual gluon fields, Eq. (3.72) then becomes 



For selfdual fields, this implies 



- (d, + igG^) 2 \*»uG$\ (1 T 75) 



(a (1 + i4 + »)^ , + = o, 



VA) = 0. 



f(+) 



has a positive spectrum we conclude that in a selfdual field i/'o, 



(3.73) 

(3.74) 

(3.75) 

(3.76) 
and thus 



and since — I d^ + igG^ 

n+ = 0. Then the index theorem (3.58) tells us that there exists exactly one left-handed (right-handed) zero 



mode per flavor in the background of an instanton (anti-instanton), and none of the opposite chirality. 



3.4 The angle 6 

In the following sections we are going to elaborate on the instanton-induced vacuum structure and the 9- 
angle which we have encountered in Section ft.l.lj . Besides the semiclassical approach and its description of 
vacuum tunneling discussed above, there are several other and partially complementary angles from which 
one can study the emergence of 9 and the implied superselection rule for non-communicating (i.e. physically 
independent) Fock spaces. Below we will consider two such approaches in more detail. Both of them lead to 
new insights into the workings of the QCD vacuum, and both of them allow the 9 structure to be derived as 
an exact, formal result, i.e. without recourse to the SCA or, in fact, any approximation. This is reassuring 
since the conditions under which the SCA works in QCD depend on the physical situation and are often 
difficult to check in the strongly coupled regime. Possibilities to learn about the consequences of the nontrivial 
gauge-group topology without relying on the SCA are therefore welcome. 

Up to now we have approached QCD in the path integral framework, which provides the most natural 
setting for the SCA. In the following we are going beyond the SCA, with the aim of augmenting our earlier 
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insights into the theta vacua 

\9) = N^e- l9n \n) (3.77) 

n 

as a familiy of (up to a phase) gauge-invariant ground states. Moreover, we will be interested in the general 
behavior of physical states under topologically nontrivial gauge transformations. It is therefore suggestive to 
chose a quantization procedure in which the Fock states play a more explicit role. Canonical quantization in 



the Schrodinger picture will be used for this purpose in Section 3.4.1 . Although instantons are indispensable 



for elucidating the tunneli ng me chanism behind the 9 structure, they will not appear explicitly in this section. 



In the subsequent Section 3.4.2 we will return to the path integral and look at the 9 vacuum from yet another 
perspective, namely that of the cluster decompositon principle. On the way, we will learn how to properly 
treat the topological charge Q of the instantons and other topologically nontrivial gauge fields in functional 
integrals. 

3.4.1 6 structure from Gaufi' law 

We begin this section with a brief summary of the Hamiltonian formulation of QCD in Minkowski space" ! 
and the canonical quantization of gauge theories. Several alternative methods for this purpose have been 
developed^], including e.g. constraint-quantization a la Dirac and the BRST formalism. In our context 
quantization in Weyl (or temporal) gauge, 

Gq = 0, Gq, ; = doGf, (3.78) 



£ = — -G^ v G a '^ v = - (d G") - -GljGlj (3.79) 



which we have already encountered in Section 3.1.1, turns out to be the most convenient. To simplify the 
discussion, we restrict ourselves to pure Yang-Mills theory (i.e. QCD without quarks) with the gauge-fixed 
Lagrangian 

iG^G Q ^ = i(d G?) 2 -J< 

since the topological effects to be discussed are rooted in the gluon sector. 
In terms of the color electric and magnetic fields 

E i = G «> B i = -^ e ^ G °jk- ( G i 3 = ~ £ i]kBl) (3.80) 

the gauge-fixed Lagrangian (in Minkowski space) is just the difference between the "kinetic" and "potential" 
gluon energy densities, 

C=\ (El - B%) . (3.81) 



The canonical variables, in Weyl gauge, are the spacial components Gf of the gauge field. In terms of those 
and the associated conjugate momenta 

" )L =G a 0i = -E? (3.82) 



1 S(doGf) 
the Hamiltonian reads 



H = J d 3 x Wd G? -£) = ±Jd 3 x (El + Bl) . (3.83) 

Note that, due to the antisymmetry of the field strength, there is no momentum conjugate to Gq, i.e. 
7Tq = SC/S (OoGq) = 0. The corresponding Yang-Mills (i.e. Euler-Lagrange) equation, Gaufi' law 

DfG b 0l = -D ab E b = 0, (3.84) 



10 A potential disadvantage of the Hamiltonian formulation is that explicit Lorentz invariance is lost because one has to single 
out a spacelike reference surface on which to impose the canonical commutation relations. 

n In QED, it can be shown that all quantization schemes lead to the same results, and one expects the same to hold for QCD. 
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is therefore just a constraint. In other words, Go does not propagate since no kinetic term for it appears 
in the Lagrangian. This is the root of several complications in the canonical quantization of gauge theories. 
One finds, in particular, that Gaufi' law, as an equation at fixed time, does not appear as one of Hamilton's 
equations. How should one then implement the full Yang-Mills dynamics in the canonical formulation? 
Certainly not by simply adding GauB' law as an operator equation: this would be inconsistent because it 
does not commute with the canonical variables. Instead, it must be imposed as a constraint on the Fock 
space, restricting the physical states to those eigenstates of H which satisfy 



DfGl |V.) = 0. 



(3.85) 



(This is analogous to how one eliminates the non-physical ghost states with negative norm, e.g., in Gupta- 
Bleuler quantization of QED or in the covariant quantization of strings.) 

In order to get a more explicit representation of the Fock space, we will now adopt the Schrodinger 
picture in which the physical states become functionals \& [G] of the canonical variables, i.e. of the spacial 
components of the gauge potential. Accordingly, the ^ [G] are eigenfunctionals of the Hamiltoniar^ 



and satisfy the Gaufi law constraint 



H^> [G] = Ef [G] 



DfG b 0i ) * [G] = 0. 



(3.87) 
(3.88) 



Let us now establish the appearance of the 9 angle as a consequence of the nontrivial topology of the 
gauge group G (i.e. ir^ (G) = Z ^ 1) in this framework |32]. The unitary gauge transformations U with 



G u 



T G„ = U (x) 



d^ + G M 



U- 1 (x) 



(3.89) 



are restricted in Weyl gauge to those "residual" ones which leave the gauge condition ( 3.78 ) intact, i.e. to 
the time-independent (and in this sense global) transformations 



u = e iu, a (a)t a = u(Lo a {x)), 



(3.90) 



where the parameters w a specify a given group element. The gauge- fixed Hamiltonian ( |3.83| ) is still invariant 
under those transformations. This implies, after quantization, that the Hamilton operator commutes with 
the operators IA (w (a?)) which furnish a unitary representation of the group on the \& [G]: 



[H,U (uj (x))] = 0. 



(3.91) 



Thus we choose the physical states ^> [G] to be simultaneously eigenstates of H and U. At first, this 
might not seem necessary since the operator Df b G\ i appearing in Gaufi' law is the generator of infinitesimal 
gauge transformationsFj, 



1 



d 3 x'w a (£') -DfG b 0i (£' ) , Gj (x) 



if) 



-Dfu b (x) = S^Gj (x) . 
t-9 



(3.93) 



Therefore, Eq. (3.88) implies that "small" gauge transformations, defined to be those which can be built 
from (up to infinitely many) infinitesimal ones and are thus continuously connected to the unit in the 



12 The canonical variables are now the classical fields C?" (x) and the corresponding momenta 

<5 



< (a?) = -Ef (x) 



SG? (x) 



(3.86) 



so that the usual canonical commutation relations are satisfied. All other (time-independent) operators, including the Hamil- 
tonian, can be constructed from the G" and 7r". 

13 Note that for small u the gauge transformation (3.89) reduces to 



G M -> G„ + — D m cj + O (u, 2 ) . 
*9 



(3.92) 
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gauge group, leave the physical states invariant. However, we have seen in Section p.l.l| that not all time- 
independent gauge transformations in QCD belong to this class. Those which do not, the so-called "large" 
(i.e. topologically nontrivial) gauge transformations, could change the phase of the physical states. 



Let us study this issue in more detail. As we have discussed in Section 3.1.1, there is in fact an enumerably 
infinite number of homotopy classes of gauge transformations (i.e. of the map S 3 — > S 3 ) which are labeled by 
an integer n. Members of the classes with n^O cannot be continuously connected to the unit element of the 
group. The classic example of a gauge transformation with n = 1 (in the group SU (2) which is embedded 
in the color gauge group SU (3)) is the "hedgehog" 



U {1) (x) = exp 



\/x 2 + . 



(3.94) 



where c is a real constant. Note that the corresponding ui?y. (x) = 27rx a /Vx 2 + c 2 do not vanish at spa- 
tial infinity (|x| — > oo), that they contain a singularity (here a branch cut), and that the limiting value 

U^ 1 ' (x) — > —1 is angle-independent. This is characteristic for topologically nontrivial gauge transforma- 
tions since the exponential can then neither be smoothly deformed to the unit element nor expanded for 
all x. A representative of the class n can be immediately obtained from LA 1 ' by using the additivity of the 
winding number under group multiplication^, 



L/ (n) (x) = \u {1) (x) 



exp 



Vx 2 + . 



(3.95) 



Thus, although all physical states are annihilated by Gaufi' law, U can have nontrivial eigenvalues for 
large gauge transformations with n ^ 0. (An explicit example is given in Appendix 0.) Since the eigenvalues 
of unitary operators must be phases, we then have 



U {n) (w (x)) * [G] = * \ u ™ MS)) g| = e i/3 " (a,(s)) * [G 



Concatenating two gauge transformations, we furthermore obtain 



U {n) {Lo)U {m) { X )^[G] 



* 



U^Uir 



G 



: U {n + m ^ (0 * [G] 



(where the group parameters lo, x, an d £ depend on x), and with (3.96) this implies 

n («) + (3 m ( X ) = (i n+m (0 ■ 



(3.96) 



(3.97) 



(3.98) 



However, the transformation behavior of ^ under gauge transformations is in fact even simpler. The existence 
of the inverse in the gauge group allows us to write the unit element as 1 = UoUq 1 for any Uq € Q with 
n = 0, so that 

U {n) {lo) * [G] = W (0) (x) U {n) (w') * [G] = e ^«(<"')tf [G] (3.99) 

i -l . 



where U^ {lo') = [U^ (%)] W'"' (to). Together with (p\96|) we then have 

(3 n (to) = f3„ (to') 



(3.100) 



and since lo' can be chosen independently of lo by selecting the appropriate Xj we conclude that the phases f3 
cannot depend on details of the g auge transformation (including its x-dependence) but only on its winding 
number. As a consequence, ( 3.98 ) simplifies to 



which finally implies 



Pn i Pm — Pn+r> 



fin = nb 



(3.101) 



(3.102) 



14 More generally, any gauge transformation U = exp [if (x 2 ) i°T ] with / (0) = and / (oo) = hi has winding number n. 
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with an arbitrary, real and state-independent angle 6 e [0, 2tt]. Altogether, we have therefore established 

U {n) («) * [G] = e m9 ^ [G] (3.103) 

for any lu and for any physical state (not just the vacuum). The angle 9, induced by the topological properties 
of the gauge group (or, more accurately, the residual, time-independent gauge transformations which preserve 
the Weyl gauge), has now appeared in a more general fashion than in our previous discussion based on the 
SCA. 



In order to link the transformation behavior ( [3.103 ) to our semiclassical tunneling picture for the 9 vacuum 



(and thus to instantons), we concentrate on the vacuum functional "Jo [G] which is characterized by the 



smallest energy eigenvalue in the physical Fock space. Furthermore, we recall from Section 3.1.1 that, while 
the vacuum gauge fields G™ before and after a gauge transformation U^ 71 ' can be continuously deformed into 
each other and are gauge-equivalent, the configurations encountered at intermediate stages of the deformation 
are generally not. If the initial and final fields correspond to the vacuum state, then at least some of the 
intermediate states must have a larger energy expectation value and thus form "potential" barriers. Those 



are penetrated by the instanton-mediated vacuum tunneling processes discussed in Section 3.1.2. And indeed, 



the resulting semiclassical vacuum state (3.77) of Bloch type realizes the above transformation law (3.103) 



U {n) \9)=N J2 z~ l8m U {n) \m) =N J2 e- l8m \n + m) = N J2 eT m ~ n) \k) = e 1 " 6 \9) (3.104) 



k=—oo 



where k = n + m. The same holds for any physical state built on \9), in accord with ( 3.103 ). 

Let us note, incidentally, that the instanton-mediated tunneling picture of the SCA provides not only in- 
sight into the mechanism which generates the ^-structure but, in suitable situations, also a convenient means 
for explicit calculations. Again, there exists an analogy with the quantum-mechanical periodic potential in, 
e.g., condensed-matter physics: the Bloch-Floquet wave function is the exact solution of the Schrodinger 
equation in the periodic potential while the tight-binding approximation provides an intuitive and effective 
tool for practical calculations. 

To summarize, we have derived the essential aspects of the QCD ^-structure, which owes its existence 
to the nontrivial topology of the gauge group SU (3), as a consequence of Gaufi' law. We have established, 
without recourse to any approximation, that QCD has a free angular parameter 9 which does not show up 
in classical chromodynamics or in the field equations and which gives rise to a superselection rule. Physical 
consequences of this new parameter will be discussed in the next section. 

3.4.2 6 via cluster decomposition 

The above discussion of the 9 angle can be complemented by following a different line of reasoning which 
starts from the cluster-decomposition requirement for QCD amplitudes 133]. As a side benefit, it gives 
new insight into the role of instantons in the path integral and in physical amplitudes. In particular, this 
approach provides an answer to the question of how the different topological charge sectors should be treated 
in functional integrals over the gluon fields. One might wonder, e.g., whether to restrict the functional 
integration to specific Q-sectors like, for example, the Q — sector in which perturbation theory takes place. 
Or do all sectors have to be included, and if so, with which weights? Interestingly, these questions can be 
answered by imposing the cluster decomposition principled which formalizes the requirement that spatially 
distant measurements yield uncorrelated results. 

At the Greens function level, the cluster decomposition principle requires that connected, time-ordered 
expectation values of products of local operators at distant positions factorize. Let us see what this implies 
for the vacuum expectation value of an operator O [q, q, G] which is composed of QCD fields and which we 
assume to be strongly localized in a spacetime volume fii. (In other words, O has support only in a small 



15 The cluster decomposition principle is one of the basic requirements for any physically sensible, local quantum field theory. 
It can be shown to hold under rather general conditions on the form of the Hamiltonian if a unique vacuum state exists [B4|| . 
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volume inside Q±.) In order to keep an open mind on the question of which Q-sectors to include, we write 

Eo- » w (Q) f VqVqDGqO [q, q, G] e -S[<i,i,GM 

where we have allowed for the most general case by summing the contributions of all topological charge 
sectors with a yet to be determined weight function w (Q) (which may be zero for some Q). Now we split the 
total Euclidean spacetime volume fl as ft = Or + f^ and note that the topological charge Q of a gluon field 
in the total volume f2 may be approximately written as a sum Q = Qi + Q2 over its topological charges Q\ 
in fii and Q2 in &,2- (Strictly speaking, Q\^ do not have to be integer since the condition ( [3.64 ) only holds 



at the boundary of f2. However, due to the localization of the topological charge density they approximately 
are.) The action is additive, too, S [0] = S [fii] + S [O2], and the functional "measure" over any field <fi on 
n factorizes as 

[■D<t>W = Y[ fd(j>{x)= \{ fd4>(xi) J[ fd<t>(x 2 )= /*P0 (Ol) x fv^. (3.106) 

Thus we have 



]T w(Q)JVGf= J2 W (Q) E ] VG q; ] E / vdSf&QAi+Q* (3-107) 

— — 00 Q — — 00 Q 1 — — oo Q2— — 00 

00 r r 

E w(Q 1 +Q 2 )JVG^ ) JVGg 2) (3.108) 



Qi,Q 

and can therefore rewrite the above matrix element as 

E Qi ,q 2 w (Ol + O2) / [VqVqDG Ql f^ O [q, q, G Ql ] e- s ^ x / [VqDqDG Q . : 1 ( ° 2) P ~ s ^ 

fain. 



M ' '° E QuQ2 ™(Qi + Q2)nVqVqVG Ql } ni e-^] x f [VqVq - VGQ2] „ e ~S[n 2 



(3.109) 

Now comes the crucial step: since the operator O is strongly localized in Cl\, cluster decomposition 
requires that (O) must be independent of what is going on in the far separated volume 0,2- This is the case 
only if 

w(Q 1 + Q 2 )=w(Q 1 )w{Q 2 ) =► w{Q) = e lQe (3.110) 

where the free angle 9 is real since w (Q) has to remain finite for all Q G [—00, 00]. As a consequence, ( |3.109| ) 
reduces to 

Eq e lQd J VqVqDG Q Q [q, q, G Q ] e ~ S ^ 

J2q etQ6 I T>qVqVG Q e- s ^ 



(0 P\ 0) n = ^ ^ J J B r^JnJ,L ■ ( 3 - m ) 



which shows that we indeed have to integrate over the gluon fields of all topological charge sectors, with a 
given, ^-dependent weight. Incidentally, the above ^-dependence is exactly what one would expect for matrix 



elements of gauge- invariant operators between the 9- vacuum states (3.77): 



T,mr,e lmf> e- m6 (m\0\n) 

(0 0==^ a - g , 1 \ (3.112) 

\ I I / J2m,n e e ( m \ n ) 

E Q ^E„(n + Q|n)] ' (6 - "> 

where Q = m — n and the expressions in the square brackets collect the matrix elements which connect 
pure-gauge sectors with winding number diffcrcnccFj Q 1 just as the functional integrals in Eq. (3.111). 



16 Recall from Eq. (3.65) that Q = n,-—,^ — Wt-^-oo is the number of potential barriers which are tunneled through. 
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Let us now interrupt for a moment the discussion of the functional integral in order to illustrate by an 
explicit example that a restriction to gluon fields of fixed Q is in conflict with cluster decomposition JJ5| . To 
this end, consider the two-point function of the pseudoscalar gluonic operator 



p(x)=G%(x)G%(x) 



(3.114) 



and restrict to the Q = sector containing an instanton and an anti-instanton of fixed size p. According to 
the cluster decomposition principle, one would expect 



n /+f (III -oo) = lim (0\Tp(x)p(0)\0) I+T =(pY 

\x\ — ►oo 



lim (0|p(s)|0) /+7 <0|p(0)|0) /+7 . (3.115) 

|cc | — >oo 



Let us check this hypothetical equation by evaluating both sides independently. For the left-hand side we 
find 



n 



1 



J+J (M-oo)= lim -[(0\Tp(x)p(0)\0) I + (0\Tp(x)p(0)\0) I ] 

\x\—^cx) L 



= - lim [(0 \p (x)\ 0) 7 (0 \p (0)| 0) 7 + (0 \p {x)\0) T (0 \p (0)| 0) f ] = A 2 

Z \x\ — >oo 

since Gj^ v = Gi tflI/ and Gj >liV = -Gj^ v imply 

(0|p(aO|0> 7 = -<0|p(aO|0> r =:A. 



On the right-hand side of Eq. ( 3.115 ) we have instead 



1 



<^=lim ~[(0\P(x)\0) I + (0\P(0)\0) I }-[(0\P(x)\0) I + (0\P(0)\0) I }=0. 



(3.116) 

(3.117) 

(3.118) 
(3.119) 



Obviously the one-instanton approximation violates cluster decomposition, and gluon fields of higher Q- 
sectors, as included in the path integral ( |3.111 ), are indispensable to restore it. 

We now continue our above line of reasoning by recasting Eq. ( [3.111 ) into a more familiar form. To this 
end, we define 

Y, VG Q = VG, S QCD = Sqcd - iQ0, (3.120) 

Q 



which implies 



(0|O|0) 



JVqVqVGO[q,q,G]e- 
J VqVqVGer S Qoi 



(3.121) 



Analyticall y con tinuing back to Minkow ski spa ce (with d A XE — id 4 XM, G a a — —iG a io, e 0123 = +1) and 
recalling Eq. ( 3.54 ), the generalized action ( 3.120 ) becomes 



D QCD,M 

which amounts to adding the term 



Sqcd, 



m 



16tt 2 



d xtr, 



(g^g^) 



Ce = -T& tr ^ G ^ u 



(3.122) 



(3.123) 



to the QCD Lagrangian (in Minkowski space). This new, renormalizable interaction had been discarded 
during the initial development of QCD since it is a total derivative (cf. (3.59)). The latter implies that it 
plays no role in perturbation theory (perturbative fields have Q — 0) or more generally for any globally trival 
gauge field, and in the field equationsPL Only the later discovery of instantons showed explicitly that this 
term can have physical consequences, the most dramatic being strong CP violation. Indeed, (3.123) breaks 

17 Hence the classical instanton solutions are unaffected by the presence of Cg, too. 



the combined charge conjugation (C) and parity (P) symmetry of the QCD L agrangian or, equivalently, its 
time-reversal (T) invariance. The latter is particularly obvious since ( 3.123 ) is at most linear in the time 
derivative. 

The physical realization of strong CP violation does not depend solely on Cg, however. This can be seen 
by invoking a flavor-dependent chiral U (1) redefinition 



If 



= ia/75 



9/ 



(3.124) 



of the quark field of flavor / in the functional integral. Due to the axial anomaly, such a transformation 
produces a nontrivial change in the measure of the quark fields, which effects a shift in the value of 9: 






(3.125) 



(This can be shown by arguments almost identical to those of Section 3.2.1.) Since the quark mass term 

Cm = ~ o Yl ["W (1 + 7s) <?/ + m* f q f (1 - 75) q f ] , (3.126) 



written here in its most general, CP and T non-conserving form with complex "mass" parameters m/ (note 
that for real m/ the 75 part vanishes and th e standard (CP-even) mass term remains), explicitly breaks 
chiral invariance, it also changes under ( 3.124 ): 



m f 



'771/ 



where we have used 



e 2M ^(l±75) = e ±2M (l±75). 



(3.127) 



(3.128) 



Since a redefinition of the path integration variables is not allowed to change physical properties, the 
latter can depend only on the invariant combination 



/ 

of 9 and the mass parameters m/, or equivalently on 



'TIki 



^arg(m/) 
/ 



(3.129) 



(3.130) 



Eq. (3.129) shows that a finite 9 would have no observable consequences (no CP violation, in particular) 
if at least one quark mass would be zero. (Although this seems quite unlikely in QCD, it cannot be firmly 
ruled out at present.) 

From the experimental bounds on the (CP-violating) dipole moment of the neutron we know that 9 has to 
be exceedingly smallr], 9 < 10~ 9 . How can the two seemingly independent terms in ( 3.130| ) cancel so (almost) 
perfectly? The unknown physical mechanism behind this cancellation and the unexplained smallness of 9 
are referred to as the "strong CP problem" . Several theoretical ideas for its solution |3(| , most prominently 
the axion models, have been proposed. Unified theories of Nature embedding QCD (e.g. string theory?) 
must even predict the value of 9 from the very first principles. This is a highly nontrivial challenge. 



8 or very close to 9 = n, which seems however to be ruled out by meson phenomenology 
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Chapter 4 

Outlook: Instantons and hadron 
physics 



We are coming to the end of these lectures and still have barely scratched the surface of instanton physics. 
This was to be expected, given the extent and variety of this field. However, in a workshop on hadronic 
physics it would hardly be appropriate to close the discussion without having given at least a glimpse of 
what impact instantons have on hadrons. 

We have argued above that instantons play an essential role in shaping the QCD vacuum. And, together 
with other intense, strongly correlated vacuum fields, they render this ground state truly complex. As a 
case in point, the "elementary excitations" of such ground states are typically not the canonical degrees of 
freedom in which we formulate the microscopic dynamics (in QCD the quarks and gluons) but rather bound 
states or "collective" degrees of freedom (in QCD the hadrons) which can be considered as disturbances of 
the vacuum "medium" . 

Condensed-matter physics supplies many interesting examples of such composite elementary excitations. 
One of them is the propagation of phonons in a crystal: the properties of phonon spectra and wave functions 
are intimately linked to the structure of the underlying crystal ground state (e.g. to its geometry, its 
distance scales, the interactions between the specific ions, etc.). As a consequence, a thorough understanding 
of excitation (phonon) properties (beyond simple mean-field approximations a la Debye) requires detailed 
knowledge of the ground state. Practically all quantum systems with many degrees of freedom share this 
requirement, and QCD is very likely no exception. Thus, knowing an important ingredient of the QCD 
vacuum wave functional, the instantons, it is natural to ask what this knowledge implies for hadron structure. 

This question is difficult to answer since interacting instanton ensembles, strongly coupled to other 
vacuum fields over large distances, do not lend themselves easily to a systematic and model-independent 
treatment (this is common to just about any infrared-sensitive problem in QCD). As a consequence, un- 
equivocal and quantitative evidence for instanton effects in hadrons turned out to be difficult to establish 
and the role of instantons in hadronic physics has remained elusive until long after their discovery. 

Nonetheless, several complementary approaches have inbetween significantly improved our theoretical 
understanding of this role. One of the first lines of attack was to include instanton-induced interactions 
into hadron models (like MIT-bag and quark-soliton models) p7| . Instanton vacuum models p, B8[ start 
at a more fundamental level and approch the physics of the instanton ensemble by approximating the field 
content of the vacuum solely as a superposition of instantons and anti-instantons. This approach, which 
neglects other (including perturbative) vacuum fields, has been developed for almost two decades and can 
describe an impressive amount of hadron phenomenology (from static properties and correlation functions 
to parton distributions f|, [3!J ) . 

As we have already noted, QCD lattice simulations recently began to complement such vacuum model 
studies by isolating instantons in equilibrated lattice configurations and by studying their size distribution 
and their impact on hadron correlators H. While results obtained from different, currently developed lattice 
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techniques have not yet reached quantitative agreement, they do confirm the overall importance of instantons 
and some bulk properties of their distribution in the vacuum. However, despite some initial attemps it will 
still take considerable time and effort before these numerically intensive simulations can establish reliable 
links to hadron structure. Moreover, lattice "measurements" usually do not give insight into the physical 
mechanisms which generate their data. For such purposes it is useful to resort to more transparent, analytical 
approaches. 

Over the last years, I have been involved in the development of such an approach, which has the capacity 
to relate the instanton component of the vacuum rather directly to hadron properties |5j . This largely model- 
independent method has led to several qualitative and quantitative insights into the role which instantons 
play in hadronic physics. Some of those I will now briefly discuss. 

The central idea is to calculate hadronic correlations functions 

ni,..., n (si, ...,x n ) = (0\TJ 1 {x 1 )J 2 (x 2 ) ■ ■ • J„(x„)|0), (4.1) 

i.e. vacuum expectation values of hadronic currents J(x) (these are composite QCD operators with hadron 
quantum numbers, e.g. Jm {x) — q(x)Tq(x) for mesons) at short, spacelike distances by means of a gen- 
eralized operator product expansion which sytematically implements instanton contributions (IOPE). The 
pivotal merit of this expansion is that it factorizes (at short distances \xi\ <C ^qcd) ^ ne contributions of all 



field modes to (4.1) into soft ones (with momenta below a given "operator renormalization scale" /u, ki < fi) 
and hard ones (ki > fi). The soft contributions, from instantons as well as from other soft vacuum fields, 
are summarily accounted for by vacuum expectation values of colorless operators (the "condensates"), while 
the hard contributions, originating from perturbative fluctuations and from small (or "direct") instantons, 
are calculated explicitly. Input are the phenomenologically known values of a few condensates and of the 



two bulk properties (3.23) and (3.24) of the instanton distribution. On this basis, the IOPE provides a 
model- independent, controlled approximation to the correlation functions at distances \x\ < 0.2 — 0.3 fm. In 
particular, it achieves a unified QCD treatment of instanton contributions in conjunction with contributions 
from long-wavelength vacuum fields and perturbative fluctuations. 

Hadron properties are obtained from the IOPE by matching it to a dual, hadronic description of the 
correlators. The latter is based on a parametrization of the corresponding spectral functions in terms of 
hadron properties and local quark- hadron duality (which is essentially a consequence of asymptotic freedom) . 
The specifics of the matching between both descriptions rely on techniques developed for QCD sum rules 
(like e.g. the use of the Borel transform) |4l| . While the application range of the IOPE approach is smaller 
than that of instanton vacuum model and lattice calculations, it has the advantages of being transparent, 
largely model-independent and fully analytical. 

Over the last years several hadronic channels have been studied in this framework, including those 
of pions [|42| [43|, baryons |4(J, 44, [45[ Q and glueballs |p7[ . In addition to quantitative predictions for 



hadron properties, these investigations have led to various qualitative insights into how instantons manifest 
themselves in hadron structure: 

1. In several hadron channels the direct instanton contributions were found to be of substantial size and 
impact. This adds to the evidence from other sources for the importance of instantons in hadron 
structure. More specifically, the results show that nonperturbative effects can strongly affect hadron 
structure already at surprisingly small distances |x| < 0.2 fm, and that most of these effects can be 
attributed to (semi-hard) instantons. 

2. Instanton effects are strongly hadron-channel selective and favor especially spin-0 meson and glueball 
channels. Moreover, various invariant amplitudes of the same correlator can receive qualitatively 
and quantitatively different instanton contributions. In such situations the neglect of hard instantons 
(which is standard practice in conventional QCD sum rules) leads to reduced stability or failure of those 
sum rules which are more strongly affected by instantons. Implementing the missing direct-instanton 
contributions led, in particular, to the resolution of long-standing stability problems in the chirally-odd 
nucleon sum- rule Q , one of the magnetic- moment sum rules p5| , the pseudoscalar pion sum rule Q 
and the (lowest moment) scalar glueball sum rule p7[ . 
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3. Stable and predictive IOPE-based sum rules exist even for correlators which do not permit a con- 
ventional QCD sum rule analysis. In particular, the first sum rule for the pion form-factor based on 
pseudoscalar currents could be established in the IOPE approach [[l2| . Its prediction for the form- factor 
agrees well with experiment in the full range of accessible momentum transfers. 

4. In the pseudoscalar pion sum rule, the direct-instanton contributions induce new, chiral-symmetry 
breaking operators which play a crucial role in generating the exceptionally small (pseudo-Goldstone) 
pion mass 143]. 

5. Instanton effects enhance the magnetic susceptibility of the quark condensate. The resulting values are 
in line with predictions of other, independent approaches. Moreover, direct instantons have a strong 
impact on one of the magnetic- moment sum rules of the nucleon and considerably improves their overall 
stability and consistency |45| . 

6. An interesting new mechanism for instanton-enhanced isospin-breaking in hadrons was found in J44| . 
Although instantons, being gluon fields, are "flavor blind" , they can strongly amplify isospin violation 
effects which originate from other, soft vacuum fields (most prominently from the difference of the 
up- and down-quark condensate) and which manifest themselves, e.g., in the proton neutron mass 
difference. Even sophisticated quark models which include instanton-induced quark interactions miss 
such effects. 

7. The scalar glueball can be (over-) bound by the instanton contributions alone J47f. In fact, this channel 
provides the first example for a sum rule which can be stabilized solely by the contributions from 
instantons. This result lends support to the findings of instanton vacuum models J48| which neglect 
the remaining contributions. Incidentally, a mainly instanton-bound ++ glueball fits naturally to the 
particularly small glueball radius tq found on the lattice J49| , which is of the order of the instanton 
size: 

r G ~ p. (4.2) 

The IOPE also provides the first set of ++ glueball sum rules which are overall consistent with the 
low-energy theorem which governs the zero- momentum limit of the corresponding correlator |J47| . 

8. Direct instantons generate by far the dominant contributions to the IOPE of the scalar glueball cor- 
relator. This makes it possible to establish approximate scaling relations between the bulk features of 
the instanton distribution and the mass mo and decay constant fa of the scalar glueball [|4j : 

m G ~ p-\ (4.3) 

So ~ ™f. (4.4) 

These relations are the first of their kind and provide a particularly direct link between instanton and 
hadron properties. 

The above findings, together with those from other sources, are beginning to assemble into a compre- 
hensive picture of how instantons manifest themselves in hadron structure and interactions. This picture 
will certainly become richer and more detailed in the future, incorporating physics ranging from nuclear 
and quark matter to hard processes, and it will very likely also teach us more about how to deal with 
nonperturbative QCD in general. 
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Appendix A 



Gauge invariance: large vs. small 



In this appendix we show that Schrodinger wave functionals which satisfy Gaufi' law may still not be invariant 
under "large" gauge transformations, i.e. those which are not continuously connected to the unit element of 
the gauge group. To this end, consider the functional 



W[G]=jd*xK» = -±- 2 j 



d xeijktr 



O, ( djG k + ^GjG k 



(A.l) 



(K^ is the "topological current" (here in Minkowski space) whose divergence we have encountered in Eq. 
(|3.59), the integrand is known as the "Chern-Simons 3-form"). One easily calculates that 



' )U fG| \e ijk tr(t a [d 3 G k + G 3 G k ]) = -Leg, (A.2) 



6Gf(x) 4tt 2 lJK v-«i"^« ' ~J~« » sir 

so that 



-D ab E b ) W [G] = -zDf ^|| = - ^ DfGfo = (A.3) 



vanishes identically as a consequence of the Bianchi identity ( |3.71 ) . This implies that W [G] varies under 



gauge transformations 5Gf = Df uj (x) with arbitrary gauge functions u> as 

where we have used Gaufi' law. Since the surface £ lies at spacial infinity we see that "small" gauge transfor- 
mations with uj a (x) — ► leave W invariant. The parameters of large gauge transformations, in contrast, do 
not vanish at spacial infinity and can therefore be topologically "active" . Under such topologically nontrivial 
gauge transformations the above surface term stays finite, 

SjjW [G] = W [ U G] -W[G]=n [U] , (A.6) 

where 

n [U] = ^J d 3 xe ljk tr [(U^diU) (U^djU) (U^dkU)] (A.7) 

provides an explicit expression for the winding number of the gauge transformation U. 
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